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1 ,  Introduction 

Air  Force  Contract  l9(f-04)-7405  was  divided  into  two  phases.  Phase 
1  was  to  review  existing  “knowledge"  in  all  aspects  of  the  atcaospheric 
environnent  and  atmospheric  physics  pertinent  to  the  problem  of  rf- 
breakdoun-.  The  results  of  this  review  (and  approximately  90%  of  the 
Final  Report)  were  published  in  book  form  and  distributed  among  the 
co-contractors  on  the  breakdown  problem.  It  was  anticipated  that  the 
review  could,  if  used  properly,  materially  increase  the  cogency  in  the 
experimental  programs  in  rf -breakdown  and  in  the  interpretation  of 
these  experiments.  Accomplishment  of  this  limited  objective  need  not 
imply  that  tf^  solution  of  the  problem  of  predicting  rf -breakdown 
thresholds  and  the  properties  of  the  resulting  plasma  for  all  types  of 
terrestrial  environments  and  at  all  altitudes  have  been  motivated.  It 
does,  however,  mean  that  the  workers  in  the  field  are  acquainted  with 
the  terrestrial  environment  as  well  as  all  the  productive  and  loss 
mechanism.s  of  electrons  -  a  state  of  awareness  which  apparently  did 
not  exist  prior  to  the  publication  of  the  review. 

Phase  2  was  an  obvicus  reaction  to  Phase  1.  In  essence,  its  objec¬ 
tive  was  to  list  ways  and  means  for  filling  some  of  the  gaps  in  existent 
knowledge  concerning  both  the  terrestrial  environment  (ambient  and 
perturbed)  and  the  man;.’  particle-particle  and  particle-field  interac - 
tiont  t  in  total  determine  the  characteristics  of  breakdown.  l!ore 
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specifically,  the  vays  and  rreans  vere  to  include  both  laboratory  (Phase 
2,  Part  A)  and  upper  attMsphere  prograss  (Phase  2,  Part  B) .  For  con¬ 
venience  of  presentation,  both  types  of  recommendations  are  included  in 
this  single  document.  The  order  or  presentation  is  selected  to  follow 
the  scheme  of  the  review  and  does  not  imply  priority  of  effort. 

As  a*result  of  the  totality  of  efforts  in  the  field  of  rf -breakdown, 
two  important  characteristics  of  the  problem  were  recognized.  First, 
the  phenomenological  theory  of  Allis,  Broxm,  etc.,  becomes  less  appro¬ 
priate  with  decreasing  pressure  -  at  least  for  those  cavity  experiments 
of  limited  size.  Moreover,  the  departure  from  theory  is  in  a  sense 
such  as  to  over-estimate  the  breakdwon  power.  For  som.e  configurations, 
it  may  be  possible  to  clarify  the  picture  by  increasing  the  dimensions 
of  the  cavity,  yet  this  has  its  inherent  drawbacks  since  there  is  a 
corresponding  increase  in  the  uncertainty  concerning  the  particular 
transmissi.n  mode  of  the  rf  perturbrtion. 

Second,  a  reliable  extrapolation  of  the  cavity  experiments  (even 
if  the  previous  objection  to  the  phenomenological  theory  can  be  elimi¬ 
nated)  must  be  considered  unlikely  in  view  of  our  relative  ignorance 
of  the  composition  of  the  upper  atmosphere  and  the  roles  of  tlie  many 
processes  such  as  solar  photo-ionization  which  can  materially  alter  the 
breakdown  characteristics.  Tlie  proper  theoretical  approach  would  be 
the  more  com.plex  Boltzmann  equation  whose  total  solution  unfortunately 
remains  beyond  the  present  state  of  the  art. 
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There  are  ways  and  means  of  improving  che  overall  posture  of  the 
state  of  knowledge  concerning  the  interactions  of  rf  energy  and  slightly 
ionized  gas  -  extending  of  course  to  the  phenomenon  of  breakdown.  It  is 
with  this  general  improvement  (in  the  state  of  the  art  In  mind)  that  a 
program  is  suggested.  The  following  list  of  items  delineates  the  broad 
prospective  of  the  program. 

(1)  Measurements  of  density,  temperature  and  composition 

(2)  Measurements  of  the  solar  flux,  absorption  and  ionization 
cross  sections 

(3)  Reaction  rates  from  ion  cloud  releases  and  from  simple 
laboratory  experiments  using  photo-ionization  and  photo¬ 
dissociation 

(4)  Photo  and  electron  impact  transition  probabilities  of 
molecular  species. 

(5)  Radiation  losses  in  ri; -excited  plasmas  obtained  by  the 
spectrometric  monitoring  of  the  Bailey,  Gordan  and  Bowles 
experiments;  and  the  theory  of  radiation  from  plasma 
inhomogeneities 

(6)  Existing  solutions  to  the  Boltzmann  equations,  the  basic 
postulates  involved,  and  the  pertinence  of  these  treat¬ 
ments  to  the  rf-breakdown  problem,  '..'here  possible, 
suggestions  are  made  for  setting  up  a  tractable  simpli¬ 
fication  of  rf -breakdown . 

(7)  A  test  of  Item  6  can  be  made  by  applying  the  theory  to 
the  Stanford  experiment  of  breakdown  in  a  spherical 
cavity  by  a  remote  rf  antenna. 
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(3)  A  direct  experir^ttaL  attack  on  the  breakdown  problem 
bj-  devising  a  roc>-et-borne  breakdown  experiirent . 

(9)  rf-breakdown  in  atmospheres  perturbed  by  nuclear  detona¬ 
tion  and  chemical  seeding. 

In  the  following  sections,  each  of  these  topics  are  discussed  in 
greater  detail . 


2 .  rp.e  Upper  Atr.osphere 


Of  the  physical  paracaters  of  the  upper  atmosphere  that  are  of 
pertinence  to  the  program  of  rf-breakdown,  the  most  important  are  density, 
composition,  and  temperature.  Each  of  these  parameters  has  been  the  sub¬ 
ject  of  considerable  study;  yet,  none  is  in  a  satisfactory  state  at  the 
present  time  -  particularly,  in  regard  to  the  requirements  of  the  rf- 
breakdown  problem. 


The  review  has  pointed  up  many  of  the  variations  in  density  and 
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the  sources  of  these  variations.  The  relation  of  density  to  such  as¬ 
pects  as  the  GSP  angle  (diurnal  variation),  the  geomagnetic  latitude 
(charged  corpuscules) ,  and  solar  activity  have  been  discussed  in  detail. 
But  most  of  these  characteristics  have  been  established  by  the  integra¬ 
tion  of  measurements  over  several  satellite  periods. 


On  the  condition  that  an  upper  atmosphere  rf -experiment  may  be 
performed  which  requires  the  local  density  (and  perhaps  the  density 
profile)  for  proper  evaluation,  this  datum  would  not  be  available. 
Furthermore,  no  reliable  data  exist  concerning  the  possible  derivations 
from  "normal  '•'■ehavior"  described  in  the  review.  Consequently,  we  advo¬ 
cate  methods  by  which  density  may  be  measured  continuously  from  the 
ground.  In  Section  2.2,  the  advantages  of  utilizing  laser  instrumenta¬ 
tion  as  the  optical  transmitter-receiver  are  given  in  detail.  .A.s  in 
th.e  case  of  the  present  optical  density  m.easuring  devices,  Ravleigh 
scatter  from  individual  molecular  and  accm.ic  species  provides  thu  basic 


m.ecnanism. 
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Tne  utility  of  the  laser  device  depends  upon  the  pulse  width  (basic¬ 
ally  the  tiir£  required  for  avalanching  of  the  photons),  the  power  trans¬ 
mitted,  and  the  receiver  sensitivity.  Pulse  width  must  be  small  if 
spatial  resolution  is  to  be  obtained,  and/or  scattering  from  the  denser 
atmosphere  eliminated  by  time -gating.  Large  transmitting  powers,  excel¬ 
lent  collimation  and  high  sensitivity  which  are  characteristic  of  the 
laser  are  exactly  the  requirements  for  a  successful  density  measurement 
via  Rayleigh  scatter. 

Above  balloon  heights,  the  composition  is  not  well  known,  and 
therefore  the  roles  of  dissociation,  mixing,  diffusion,  etc.,  have  been 
the  subject  of  considerable  conjecture.  The  problem  is  inherently  more 
difficult  than  that  of  the  density,  temperature,  etc.,  particularly 
since  trace  species  may  be  of  considerable  importance.  If  laser  instru- 
mentation  can  be  built  with  sufficient  power,  sensitivity,  and  resolu¬ 
tion,  it  may  be  possible  to  use  this  device  in  composition  studies. 

The  basic  idea  would  be  to  supplement  ordinary  non-resonant  scatter 
with  resonant  scatter  by  the  appropriate  choice  of  the  frequency.  This 
m.a'/  mean  selecting  frequencies  in  the  infrared  or  near  UV,  and  possibly 
a  fundamental  extension  of  the  principle  of  stimulated  em.ission.  The 
background  may  also  limit  the  tool  to  night-tim.e  use. 

A  more  immediately  available  instrument  is  the  mass  spectrometer 

(herzog  version  of  the  Redhead  Gauge)  presently  being  used  to  get  the 

He  density  profile.  The  instrument  can  be  used  equally  well  for  any 
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specie  and  has  an  absolute  lim.it  of  10  m.olecules  cm.  ,  or  at  atmos- 
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pheric  pressure  can  detect  a  trace  material  »„hosa  density  is  10  of 
the  total.  Of  course,  the  instrument  must  be  rocket  or  satellite -borne ; 
this  is  perhaps  less  convenient  and  more  expensive  than  ground-based 
measurements.  Details  of  the  system  are  given  in  Section  22- 

2,1  Lasers  for  Continuous  Density,  Composition  and  Temperature 
Measurements* _ 


In  the  past,  the  scattered  light  from  searchlight  beams  has  been 
used  in  measurements  of  upper  atmospheric  density.  The  recent  develop¬ 
ment  of  lasers  and  masers  (sources  of  highly  monochromatic  coherent 
radiation  of  high  total  power  and  extremely  high  power  per  cycle  per 
steradian)  has  made  feasible  their  utilization  as  geophysical  probes. 
Lasers  may  be  employed  in  making  the  following  measurements: 

(a)  Measurements  of  density  and  hence  temperature  as  a  func¬ 
tion  of  height  to  75  fcms  and*possibly  higher  to  100  kms. 

(b)  Measurements  of  height  distribution  of  various  atmospheric 
species  via  e.g.,  a  laser  in  the  region  2500-3000S  for 
ozone  or  an  iraser  in  an  appropriate  IR  region  for  other 
species  (CO^) . 

(c)  Possible  measurement  of  winds  at  extreme  altitudes  by 
means  of  the  Doppler  shift  of  the  extremely  narrow  wave- 
bond  . 

(d)  Possible  mearure-ren:  of  temperature  through  Doppler  shift. 

(e)  Measurement  of  dust  layers  and  micro-meteorites , 

’■Taken  in  part  from.  J.  Pressman  (aCA)  . 
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2.1  1  Description 

The  laser  as  a  source  has  several  unique  aspects.  First  of  all, 
light  emitted  is  coherent  and  this  coherence  is  both  latitudinal  and 
longitudinal.  So  then  the  laser  acts  In  the  manner  similar  to  a  radio 
transmitter.  This  coherence  arises  fundamentally  out  of  the  nature  of 
the  induced  emission  process.  A  second  major  property  is  that  of  mono¬ 
chromaticity;  this  monochromaticity  is  such  that  the  line  emitted  has  a 
width  much  less  than  the  natural  width.  A  third  important  feature  is 
the  very  narrow  angle  of  divergence  of  the  beam.  The  final  feature  is 
the  large  amounts  of  light  power  that  can  be  generated. 

Amplifiers  and  oscillators  using  atomic  and  molecular  processes,  as 
do  the  various  varieties  of  masers,  may  extend  far  beyond  the  frequency 
range  which  has  been  generated  electronically  into  the  infrared,  the 
optical  region,  or  beyond.  Such  techniques  have  realized  the  attractive 
promise  of  coherent  amplification  at  these  high  frequencies  and  of  gen¬ 
eration  of  very  monochromatic  radiation.  As  it  was  attempted  to  extend 
maser  operation  towards  very  short  wavelengths,  a  number  of  nexj  aspects 
and  problems  arose,  which  required  quantitative  reorientation  of  theory 
and  considerable  modification  of  che  experimental  techniques  used. 
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the  fol Icv.-ir.g  discussion,  roughly  reasonable  values  of  design 
rs  will  be  used.  For  comparison,  the  characteristics  consider^ 
aper  of  Schawlow  and  Townes''”  ''  will  be  presented  for  masers 
g  in  the  normal  microwave  range.  Here  an  unstable  ensemble  of 
r  molecular  system.s  is  intrcduced  into  a  cavity  v;hich  has  one 
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resonant  node  near  the  radiative  transitions  of  these  systems.  Such  an 
ensemble  may  be  located  in  a  wave  guide  rather  than  in  a  cavity  but 
again  there  would  be  typically  few  modes  of  propagation  allowed  by  the 
wave  guide  in  the  frequency  range  of  interest.  The  condition  of  oscilla¬ 
tion  of  n  atomic  systems  excited  with  random  phase  and  located  in  a  cavity 
of  appropriate  frequency  may  be  written  • 

n  >  hV  A  v/(4m  Q^)  (2-1) 

where  n  is  more  precisely  the  difference  Nj^-n^  in  the  number  of  systems 
in  the  upper  and  lower  states,  V  is  the  volume  of  the  cavity,  A  v  is  the 
half-width  of  the  atomic  resonance  at  half-maximum  intensity,  assuming  a 
Lorentzian  line  shape,  p,  is  the  matrix  element  involved  in  the  transition, 
and  is  the  quality  factor  of  the  cavity.  ^ 

The  energy  emitted  by  such  a  maser  oscillator  is  extremely  mono¬ 
chromatic,  since  the  energy  produced  by  stimulated  emission  is  very  much 
larger  than  that  due  to  spontaneous  emission  or  the  normal  background  of 
thermal  radiation.  The  frequency  range  over  which  appreciable  energy  is 
distributed  is  ^innroxifT’ately  by 

5v  -  4,-  kT  (A  v)“  /  P  (:2-2) 

where  A  v  is  the  half-width  at  ha  If -maximum  of  the  resonant  response  of 
a  single  atomic  system,  P  is  the  total  power  emitted,  k  is  Boiczm.ann's 
constant,  and  T  the  absolute  temperature  of  the  cavity  walls  and  wave 
guide  Since  in  all  maser  oscillators  at  microwave  frequencies. 
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P  >  >  kT  i  1/,  the  radiation  is  largely  emitted  over  a  region  very  much 
smaller  than  A  Vt  or  5v  <  <  A  v- 


The  necessary  condition  for  oscillation  nay  be  obtained  by  requiring 
that  the  power  produced  by  stimulated  emission  is  at  least  as  great  as 
that  lost  due  to  all  loss  processes.  That  is. 


u'E  ^  h  V  n  >  V 
n  4n  A  i/  *”  8jt  t 


(2-3) 


2 

where  is  the  matrix  element  for  the  emissive  transition,  E  is  the 

2 

mean  square  of  the  electric  field  (for  a  multi-resonant  cavity,  E  may 
be  considered  identical  in  all  parts  of  the  cavity),  n  is  the  excess 
number  of  atoms  in  the  upper  state  over  those  in  the  lower  state,  V  is 
the  volume  of  the  cavity,  t  is  the  time  constant  for  the  rate  of  decay 
of  the  energy,  and  A  v  is  the  half-width  of  the  resonance  at  half  maximum 
intensity,  if  a  Lorentzian  shape  is  assumed.  The  decay  time  t  may  be 
written  as  2:r  v/Q,  but  may  also  be  expressed  in  terms  of  the  reflection 
coefficient  a  of  the  cavity  walls. 


t  -  6V  /  (l-a)  Ac 


(2-4) 


where  A  is  the  wall  area  and  c  the  velocity  of  light.  -For  a  cube  of 
dimension  I,  t  -  L  /  (l-CX)  c.  The  condition  far  oscillation  from 
Eq.  (2-3;  is  then 


n  ^ 


h(l-a) 
16  ^  4 


Ac 

2 


(2-5) 


At  low  pressure,  most  infrared  or  optical  transitions  will  have  a 
width  A  V  determined  fay  Doppler  affects.  The.n  the  resonance  half-width 
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where  m  is  the  molecular  mass,  k  is  Boltzmann's  constant,  and  T  the 
temperature.  Because  of  the  Gaussian  line  shape  in  this  case,  Eq,  (2-5) 


becomes 


>  A  V  h  (l-Cit)Ac 

16ic^  (itln2)^ 


>  h(l-a)A  /  2W 

"  “  2  2  1  Jtm 

16n  p.  » 


(2-7) 


(2-8) 


It  may  be  noted  that  Eq.  (2-8)  for  the  number  of  excited  systems 

required  for  oscillation  is  Independent  of  the  frequency,  Also,  this 

number  n  is  not  impractically  large.  If  the  cavity  is  a  cube  of  1  cm 

•»  18  o 

dimension  and  a  -  0.98,  p  ■:!  5x10  esu,  T  -  400  K,  and  m  100  amu, 
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one  obtains  n  =  5x10  . 


The  minimum  power  which  must  be  supplied  in  order  to  maintain  n 
systems  in  excited  states  is 


P  T  nhv/x 


This  expression  is  independent  of  the  lifetime  of  Che  excited  species. 
However,  one  must  consider  that  if  there  are  .alternate  modes  of  decay  of 


each  system,  as  by  collisions  or  other  transitix 


the  necessary  power 


may  be  larger  than  that  given  by  Eq.  i.2-9)  an.l  dependent  on  details  cf 
the  system  involved.  Furthermore,  some  quantum  of  higher  frequency  than 
that  emitted  will  normally  be  required  to  excite  the  system,  which  xvill 
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increase  the  power  somewhat  above  the  value  given  by  Eq.  (2-9) .  Assuming 

the  case  considered  above,  i.e.,  a  cube  of  1  cm  dimension  with  a  -  0.98, 

4 

A  -  10  A,  and  broadening  due  to  Doppler  effect,  Eq.  (2-9)  gives 
-3 

P  ^  0.8x10  watt.  Supply  of  this  much  power  in  a  spectral  line  does 
not  seem  to  be  extremely  difficult  and  obviously  has  already  been 
accomplished. 

The  power  generated  in  the  coherent  oscillation  of  the  maser  may  be 
extremely  small,  if  the  condition  of  instability  is  fulfilled  minimally, 
and  hence  can  be  much  less  than  the  total  power,  which  would  be  the  order 

-3 

of  10  watt,  radiated  spontaneously.  However,  if  the  number  of  excited 
systems  exceeds  the  critical  number  appreciably,  then  the  power  of  stimu¬ 
lated  radiation  is  given  roughly  by  hv  times  the  rate  at  which  excited 
systems  are  supplied,  if  the  excitation  is  not  lost  by  some  othef'*processes , 
The  electromagnetic  field  then,  builds  up  so  that  the  stimulated  emission  may 
be  appreciably  greater  than  the  total  spontaneous  emission.  For  values  even 
slightly  above  the  critical  number,  the  stimulated  power  is  of  the  order  of 
the  power  nhv/t  supplied,  or  hence  of  the  order  of  one  milliwatt  under  the 
conditions  assumed  above. 

The  most  obvious  technique  for  supplying  excited  atoms  is  excitation 
at  a  higher  frequency,  as  in  optical  pumping  or  a  three-level  maser  system 
The  power  supplied  must  be  larger  than  the  emitted  power  in  Eq.  (2-8).  It 
is  not  necessary  that  the  pumping  frequency  be  much  higher  than  the  fre¬ 
quency  emitted,  as  long  as  the  difference  in  frequency  is  much  greater  than 
kT/h,  which  can  cissure  the  possibility  of  negative  temperatures"  Since,  for 
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the  high  frequencies  required,  an  incoherent  source  of  puraping  power 
must  be  used,  a  desirable  operating  frequency  would  be  near  the  point 
where  the  most  quanta  are  emitted  by  a  given  transition  in  the  discharge 
or  some  other  source  of  high  effective  temperature  This  maximum  will 
occur  somewhere  near  the  maximum  of  the  blackbody  radiation  at  the  effec¬ 
tive  temperature  of  such  a  source,  and  most  often  in  the  visible  or 
ultraviolet  region.  The  number  of  quanta  required  per  second  would 
"  probably  be  about  one  order  of  magnitude  greater  than  the  number  emitted 
at  the  oscillating  frequency,"  so  that  the  input  power  required  would  be 
about  ten  times  the  output  given  by  Eq,  (2-9),  or  ten  milliwatts.  This 
amount  of  energy  in  an  individual  spectroscopic  line  is,  fortunately, 
obtainable  in  electrical  discharges  and  this  technique  has  been  developed 
by  Javan . 

Monochromaticity  of  a  maser  is  very  intimately  connected  with  its 
amplifier  noise  properties.  If  there  is  considered  first  a  maser  cavity 
for  optical  or  infrared  frequencies  which  supports  a  single  isolated  mode, 
then  as  in  the  microwave  case,  it  is  capable  of  detecting  one  or  a  few 
quanta,  corresponding  to  a  noise  temperature  of  hu/k.  However,  at  a 
wavelength  of  one  micron,  this  noise  temperature  is  about  14,000*^K,  and 
hence  not  remarkabl/  low.  Furthermore,  other  well-known  quantum  detectors, 
such  as  a  photoelectric  Cube,  are  capable  of  ;!etecting  a  single  quancum- 
At  such  frequencies,  a  maser  has  no  great  advantage  over  well-known 
techniques  in  detecting  small  numbers  of  quanta. 

There  is  now  examined  Che  extent  to  which  the  normal  line  width  of  the 
emission  spectrum  of  an  atomic  system  will  be  narrowed  by  maser  action,  or 
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hence  how  monochrcmatic  the  emission  frsa  an  infrared  or  optical  maser 
would  be.  Analysis  was  made  by  Schawlow  and  Townes  of  the  number  of 
excited  Systems  required  to  produce  stimulated  power  which  would  be  as 
large  as  spontaneous  emission  due  to  that  of  a  multimode  cavity  whose 
frequencies  lie  within  the  resonance  "width  of  the  system.  They  assumed 
for  the  moment  that  a  single  mode  can  be  isolated.  Spontaneous  emission 
into  this  mode  adds  waves  of  random  phase  to  the  electromagnetic  oscilla* 
tions,  and  hence  produces  a  finite  frequency  width  which  may  be  obtained 
by  analogy  with  Eq.  (2-2)  as 

i!i  V  (4jt  hv/P)  (A  w)^  (2-10) 

osc 

where  A  v  is  the  half-width  of  the  resonance  at  half-maximum  intensity, 
and  P  the  power  in  the  oscillating  field.  Note  that  kT,  the  energy  due 
to  thermal  agitation,  has  been  replaced  fay  hv,  the  energy  in  one  quantum. 
Usually  at  these  high  frequencies,  hv  »  kT,  and  there  is  essentially  no 
"thermal"  noise.  There  remains,  however,  "zero-point  fluctuations"  which 
produce  random  noise  through  spontaneous  emission,  or  an  effective  tempera¬ 
ture  of  hv/k. 

For  Che  case  considered  nunierlcally  above,  4:t  hv  A  '.  /P  is  near  10  ^ 

when  P  is  given  by  Eq ,  (2-9),  so  that  L.  v  ~  10  ^  A  v ,  This  corresponds 

osc 

CO  a  remarkably  monochremac ic  emission  However,  for  a  multimode  cavity, 
this  very  monochromatic  emission  is  superimposed  on  a  background  of  stimu¬ 
lated  em.isaion  which  has  w'id_th  a  ,  and  which,  for  the  power  P  assumed,  is 
of  intensity  equal  to  that  of  the  stimulated  emission  Only  if  the  power 
is  increased  by  same  additional  ractar  cf  about  ten,  or  if  the  desired  mode 
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is  separated  from  the  large  number  of  undesired  ones,  would  the  rather 
monochromatic  radiation  stand  out  clearly  against  the  much  wider  fre¬ 
quency  distribution  of  spontaneous  emission. 

The  above  analysis  has  been  given  in  some  length  to  indicate  the 
technique  of  construction  of  the  laser  so  as  to  have  an  adequate  concept 
of  the  power  range,  wavelength  range  and  degree  of  purity  of  the  light 
emitted.  This  analysis  follows  closely  Schawlow  anti  Townes*  development 
which  is  basic  to  an  understanding  of  the  physics  of  laser  operation. 

2.1.2  Summary  of  Existing  Laser  "Hardware"  and  Potential  Developments 

The  preceding  paragraphs  have  delineated  the  physics  of  the  laser 
device  as  part  of  the  necessary  background  for  the  uses  to  be  discussed 

•• 

later.  It  is  the  function  of  this  section  to  present  a  broad  view  of  the 
present  .state  of  art  in  this  field  as  well  as  its  potential  development 
so  33  to  place  ourselves  in  an  adequate  position  to  discuss  the  design 
of  laser  geophysical  probe  devices. 

2.1  2,1  Laser  "Transmitters" 

First, as  a  preliminary  set  of  remarks,  it  must  be  said  that  the  laser 
field  is  moving  forward  e.xtremely  rapidly.  Steady  improvements  are  being 
mad-e  in  power  emitted,  range  of  wavelengths  available,  modulation  capability, 
and  total  system  efficiency  At  present,  one-half  dozen  lasers  are  available 
as  aff-the-sheir  items.  In  Table  2-1  are  given,  for  example,  the  character¬ 
istics  of  "Vireo"  lasers,  presently  commercially  available.  In  addition, 
maty  laser  components  such,  as  crystals,  Fabry-Perot  end-plate,  light  sources 
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for  optical  pumping,  etc.,  are  commercially  available  for  use  in  the 
construction  of  special  laser  devices  designed  to  nseec  specific  needs. 

Present  ruby  masers  have  produced  pulsed  optical  beams  of  10  kilo¬ 
watts  power  or  still  shorter  pulses  of  perhaps  100  kilowatts  peak  power. 
This  radiation  is  concentrated  in  a  band-width  of  about  0.02  cm  ^ . 

Another  type  of  maser,  operating  continuously  at  about  0.02  watts,  emits 
a  wave  which  has  been  shown  to  be  in  phase  over  the  entire  maser  reflector 
surface  and  to  be  concentrated  in  a  frequency  interval  of  about  10  kilo¬ 
cycles.  There  seems  to  be  no  general  reason,  other  than  the  necessary 
dissipation  of  power,  why  solid  state  optical  masers  cannot  operate  con¬ 
tinuously  at  high  power  and  with  a  short-time  monochromaticity  close  to 
theoretical  expectation  or  hence  with  a  frequency  width  very  much  less 
than  1  megacvcle/sec . 

If  a  laser  produces  a  wave  of  wavelength  X  with  constant  phase  over 
a  surface  of  diameter  d,  the  angular  width  of  the  radiating  beam  is 
approximately  X/d.  This  can  be  reduced  still  further  by  use  of  an  auxil¬ 
iary  optical  system  such  chat  the  angular  width  may  be  reduced  to  X/D 
where  D  is  equal  Co  the  diameter  and  focal  length  of  the  lens  The  con¬ 
clusion  here  is  that  the  angular  divergence  of  the  beam  is  limited  only 
by  th.e  optical  distortions  of  the  system  and  not  by  the  nature  of  the 
source 

For  example,  the  above  10  kilov.att  continuous  source  would  give  with 
a  200"  reflector  a  beam  of  intensity  (according  to  Townes) 
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I  ^  Flux/(X/D)  =  10  watts/starad  at  5000  A,  or  per  wave  number 

22  -1  '5  *1 

I  -  3x10  watts/sterad  cm  for  a  bandwidth  of  A  v  ^  3x10  cm 

V 

In  addition  to  the  chromium  doped  ruby,  lasers  are  also  being  developed 
using  samarium  and  uranium  doped  calcium  fluoride.  These  give  a  variation 
in  wavelength  and  are  more  adaptable  to  CW  operation.  At  present,  solid 
state  continuous  wave  operation  is  on  the  verge  of  actuality. 

Continuous  wave  operation  by  means  of  gas  discharge  has  already  been 
f2-2) 

achieved  by  Javan  at  Bell  Telephone  Labs  albeit  at  relatively  low 

energies .  At  present,  systems  are  already  being  designed  in  the  megawatt 
to  10  megawatt  output  region  by  various  groups  for  pulsed  operation. 

2. 1.2. 2  Laser  Receivers 

As  mentioned  earlier,  some  conventional  devices  such  as  photomulti¬ 
pliers  have  the  capability  of  detecting  a  few  quanta  so  that  no  great 
advantage  is  gained  here  for  the  use  of  a  laser  receiver.  However,  the 
major  advantage  of  such  a  receiver  is  that  it  rejects  to  an  extremely 
large  degree  extraneous  background  light  and  also  provides  a  nearly 
noiseless  receiver.  This  would  hold  true  even  when  such  a  receiver 
would  be  aligned  v;ith  the  sun. 

Such  receivers  are  not  presently  on  the  market  and  commercially 
available.  However,  such  receivers,  it  is  expected,  will  be  announced 
in  the  next  few  months.  This  instrument  would  be  very  sim.ilar  ir  con¬ 
struction  to  that  of  the  transmitterT  In  the  receiver,  in  contradistinction 
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to  the  transnitter,  the  Fabry-Perot  end  parallel  mirrors  would  be 
slightly  transparent.  The  major  problem  here  is  that  the  receiver  must 
be  exactly  matched  to  the  transmitter,  i.e. ,  the  receiver  mirror  spacing 
must  be  in  resonance  for  the  transmitter  signal.  Additionally,  the  axis 
of  both  receiver  and  transmitter  must  be  parallel  within  a  very  close 
tolerance  to  reject  extraneous  background  light. 

The  internal  noise  of  the  receiver,  due  to  its  owi  omni-directional 
glow,  is  considered  negligible.  Also,  thermal  noise  is  considered  insig¬ 
nificant  because  the  receiver  operates  at  a  temperature  below  that 
necessary  to  emit  energy  in  the  visible  portion  of  the  spectrum.  Hence, 
the  laser  receiver  has  a  very  definite  advantage  over,  say,  a  photo¬ 
multiplier  in  this  respect.  Such  an  advantage  may  be  minimized  if 
photomultipliers  are  used  in  conjunction  with  very  narrow  band  filters. 
But  such  accessory  filters  cannot  be  expected  to  match  the  built-in 
selective  natural  monochromatic  response  of  the  laser  receiver, 

2.1.3  Scattering  of  a  Searchlight  Beam 

Because  the  Rayleigh  scatter  of  a  searchlight  beam  offers  a  simple 
but  effective  method  of  measuring  the  density  and  temperature  of  the 
upper  atmosphere,  it  would  be  desirable  if  such  a  method  could  be 
extended  to  higher  altitudes  (above  60  kra)  by  means  of  a  laser  beam. 

The  elements  of  the  scattering  technique  are  presented  in  the  following 
sections 

The  light  that  is  reflected  by  the  upper  atmosphere  into  the 
receiver  is  proportional  to  the  number  of  scattering  centers  present 


at  the  altitude  of  interest.  It  is  assumed,  in  the  discussion  that 
follows,  that  the  composition  of  the  atmosphere  remains  constant.  Pre¬ 
vious  results  confirm  that  the  method  of  scattering  can  be  described  by 
the  Rayleigh  law,  in  which  it  is  assumed  that  the  particles  have  diameters 
that  are  small  compared  to  the  wavelength  of  the  light. 

Johnson  and  Jones  have  estimated  the  shape  and  magnitude  of  the 
reflected  signal.  The  discussion  that  follows,  which  obtains  their 
result,  leans  heavily  on  their  presentation. 

Data  on  the  upper  atmosphere  are  available  from  the  Rocket  Panel, 
0-3) 

Elterman,  and  others.  The  magnitude  of  the  pressure  and  ten^ra- 

ture,  up  to  at  least  100  km,  can  be  deduced  from  their  results.  Bence, 
it  is  possible  to  calculate,  using  the  standard  gas  laws,  the  number  of 
molecules  per  cc  as  a  function  of  altitude. 


Hulbert  has  given  a  convenient  expression  for  the  amount  of  light, 
in  lumens,  that  is  scattered  back  per  unit  solid  angle,  per- cm  length  of 
path,  at  an  angle  of  5  to  the  incident  beam. 


i 
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io  2n-  (^^-l)^(Ucos“o) 


(2-11) 


where  X  is  the  wavelength  of  light,  the  refractive  index  at  wavelength 
\,  n  the  number  or  molecules  per  cc  and  i^  the. incident  beam  flux  in.  lumens 


We  also  have 
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with  a  value  of  1.08x10  for  OL  as  given  by  Johnson,  Meyer,  Hopkins 
(2-4) 

and  Mock.  Since  the  variation  of  (X^  with  wavelength  is  small,  its 

value  for  these  calculations  may  be  taken  to  be  a  constant. 

The  Rayleigh  formula  is  known  to  be  incorrect  by  a  factor  of  two. 

If  all  of  the  above  relations  are  combined,  we  find  that  the  intensity 
of  the  back  scattered  light,  in  lumens,  for  O-O,  per  lumen  incident,  per 
cm  length  of  beam,  per  unit  solid  angle,  is 

i  =  9.05  X  10"^®  X  n  (2-13) 

s 

The  Intensity  of  light  reaching  any  particular  altitude  is  equal  to 
the  intensity  of  the  light  source  multiplied  by  the  transmission  factor 
up  to  the  altitude.  Similarly,  the  intensity  of  the  back  scattered  light 
reaching  the  receiver,  is  equal  to  the  intensity  of  the  light  scattered 
in  the  backward  direction  multiplied  by  the  same  transmission  factor. 


The  transmission  losses  are  due  to  two  factors.  The  first  is  due 
to  molecular  scattering,  and  applies  for  all  altitudes  The  second  is 
due  to  large  particle  scattering  and  applies  only  to  the  low  altitudes. 

The  transmission  losses  that  are  due  to  molecular  scattering  can  be 

computed  from  the  Rayleigh  formula 

3  ^ 

32n  a  n 


3.\ 


(2-14) 


which  is  obtained  from  Eq.  (2-11)  by  integration  over  a  sphere.  Upon  so 
doing,  Johnson  and  Jones  conclude,  taking  into  account  losses  in  the 
lover  atmosphere,  that  the  transmission  losses  are  not  great  and  result 


21 


in  a  total  transmission  of  about  85%  for  the  first  50  km  while  the 
losses  above  50  km  are  negligible. 

The  intensity  of  the  back  scattered  light  entering  the  receiver 
searchlight  is  given  by 

i^  =  const  n  x  A/h^  (2-15) 

where  A  is  the  area  of  the  receiving  mirror  and  h  is  the  altitude. 

Utilizing  the  availabl^data  for  n,  the  size  of  the  searchlight,  etc.,' 

a  plot  of  the  expected  scattered  signal  may  be  obtained.  Figure  2-1, 

(2-5) 

taken  from  the  paper  of  Johnson  and  Jones,  shows  the  estimated 

signal. 

If  negligible  losses  in  the  searchlight  beam  above  10  km  are 
assumed,  then  the  ratio  of  the  signal  for  any  two  altitudes  is 

n  n^  (S/S^)  (h/h^)^  (2-16) 

where  S  represents  the- magnitude  of  the  signal  from  altitude  h  at  the 
receiver  If  the  lover  altitude  is  accessible  to  the  radiosonde,  it 
then  becomes  possible  to  compute  the  density  profile,  with  the  radiosonde 
altitude  as  the  reference  altitude. 


It  Is  important  to  note  for  further  reference,  chat  between  15  and 
60  km  there  is  a  change  in  the  light  signal  covering  four  to  five  orders 
of  magnitude;  whereas,  from  the  ground  to  60  km  the  change  in  the  light 
signal  covers  10  orders  of  magnitude  as  can  be  seen  from  Fig.  2-1  In 


directly  the  P.avleigh  scatter  "reflectivitv"  for  ba 
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Figure  2-2.  Atmospheric  Reflectivity 
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scatter  as  a  function  of  altitude  and  in  Table  2-2  the  data  for  4000 
and  2000  A  are  tabulated. 

In  Fig.  2-3  is  given  an  example  of  Elterman’s  results  for  density 
and  temperature. 

2.1.4  Analysis  of  Laser  Possibilities  as  Geophysical  Probe 

The  purpose  of  this  section  is  to  indicate  the  geophysical  poten¬ 
tialities  of  the  use  of  lasers  in  order-of-magnitude  fashion.  First, 
there  will  be  indicated  its  use  in  measuring  densities. 

2. 1.4.1  Measurement  of  Densities 

To  obtain  an  indication. of  the  potential  use  of  lasers  for  this 
purpose  a  simple  line  of  argument  will  be  used.  Elterman  in  his  experi¬ 
ments  used  as  a  source  a  carbon-arc  lamp.  The  estimated  total  flux  in 
radiant  energy  was  about  100  watts  spread  over  about  6000  A-  The  flux 
for  lasers  continuously  operated  can  run  to  about  10  kilowatts.  This 
gives  a  factor  of  increase  of  about  ICO  in  light  output.  The  normal 
atmospheric  density  decreases  exponentially  with  a  relaxation  height  of 
approximately  8  kms  for  a  factor  of  e.  This  gives  an  additional  height 
for  probing  of  about  36  kilometers;  so  that  since  Flterman's  experiments 
routinely  gave  data  up  to  60  kms,  the  lasers  of  the  type  mentioned  should 
go  up  to  som.e  36  km.s  additional  or  to  about  95  kms. 

f 

Furthermore,  the  carbon  arc  is  broad-band,  about  6000  A  vhile  the 
laser  has  a  width  less  than  1  A.  This  reduction  in  bandwidth  allcws 
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AT>JOSPHERIC  REFLECTIVITY  (RAYLEIGH) 


Altitude 
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kiloneters 

N  cm 

Tj  cm  (0.4^) 

TJ  cm  (O.Tji) 

300' 

1.52  X  10® 

-19 

2.87  X  10 

4.59  X  10‘^® 

200 

2.95  X  10^ 

5'.  57  X  lO'^® 

8.92  X  10‘^^ 

• 

100 

13 

1.49  X  10 

-14 

2.82  X  10  ^ 

-13 

4.50  X  10 

- 

• 

50 

2.23  X  10^® 

4.22  X  lO"^^ 

6.74  X  10"^° 

• 

25 

17 

8.32  X  10 

1.57x10'^ 

2.53  X  lO"® 

11 

7.57  X  10^® 

1.43  X  lO'® 

2,29  X  lO"^ 

5 

19 

1.54  X  10 

2.91  X  10*® 

4.65  X  lO"^ 

0 

2.55  X  10^^ 

4.82  X  lo'® 

7.70  X  lO'^ 

use  of  a  narrowband  -filter  or  laser  receiver  which  should  cut  down  the 
background  illumination  from  the  airglow  by  a  factor  of  1000.  The  exceed¬ 
ingly  narrow  cone  of  divergence  of  the  laser  which  is  less  than  15’  of  arc 
gives  when  compared  to  Elterman’s  instrument,  which  had  a  divergence  of 
1.5  ,  a  factor  of  geometrical  narrowing  of  approxirsately  100,  This  should 
also  improve  the  background  noise  by  a  large  factor. 

The  improvement  in  background  by  this  geometrical  narrowing  is  due 
to  two  main  factors.  First,  the  natural  background  due  to  airglow. 
aodiacal  light,  star  light,  etc.  is  diminished  directly  by  this  narrowing, 
secondly,  because  of  the  smaller  acceptance  cone  the  effect  of  secondary 
scattering  from  the  laser  light  itself  is  diminished.  In  the  CW  case  this 
is  a  direct  effect.  In  the  pulsed  case  this  effect  is  due  to  a  diminution 

of  taii-off  scattering  which  obscures  the  returning  pulse  from  the  high 
altitudes. 

Also,  because  of  the  better  height  resolution  due  to  the  smaller 
volume  anticipated  by  the  transmitter  and  receptor  cones,  the  measure¬ 
ments  of  temperatures  through  the  hydrostatic  equation  can  be  more  e.xact. 

As  a  check  on  the  analogous  argument  via  Elter=:an’s  experiment 
another  computation  is  made  ab  initio.  We  have  ' 


F,  =  total  flux  at  receiver 

~  initial  flux  from  transmitter 


Ti  -  Rayleigh  backscatter  from  a  one  cm^  volume 
=  length  of  scaccering  volume 
A  =  cross  section  of  receiver 
^  ~  height  of  scattering  volume 

We  chose  a  setup  similar  to  Eiterman's.  Let 

h  =  100  fcms  =  10^  cms  *  10  kilowatts  at  3000  A  =.  2x10^^  photons/sec 
1  kffl  -  10  cms,  A  —  3xi.O  cm  (area  of  a  two-meter  receiving  mirror) 

and 

^100  ~ 

Hence, 

F  _  j-XlO  ^^Xl0^x3xl0^x2xl0^^  4 

^  ^q14  ~  photons  -  a  detectable  amount  . 

The  optimal  technique  to  be  used  should  be  investigated.  Some  of 
the  systems  are 


1-  Continuous  emission  with  crossed  beams 

a.  steady  signal 

b.  scatter  modulaoei 

c.  rotating  beam  -oduLated 
2  Pulsed  emission 

o-  total  count  of  return 
b.  envelope  analysis 


Ke.nce,  the  seccnd-generat: 
scattering  up  to  conserva: 
iCO  kms 


■2-iT5  should  measure  density  via  Rayleigh 
73  kns,  and  probably  higher  up  to  possibly 
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A  similar  argument  holds  for  pulsed  lasers  which  can  run  up  to 

-3  ■  , 

10,000  megalumens  (lumen  =  1.61x10  watts  at  X  of  maximum  visibility) 

/2«g\ 

compared  to  the  50  megalumens  previously  used  by  Friedland  et_al.  ' 
Also,  the  pulse  length  can  be  much  shorter  and  hence  minimizes  the  tail- 
off  problem  that  Friedland  ran  into.  Hence,  substantial  range  in  altitude 
probing  seems  probable  over  his  results. 

2. 1.4. 2  Measurement  of  Temperature 

By  use  of  the  hydrostatic  equation  it  is  possible  to  compute  the 
temperature  from  the  density  profile  in  a  manner  similar  to  previous 
attempts.  The  temperature,  as  a  function  of  density,  is: 

T  =  ^/R  p  P  gpdh- 

WW 

where  p  «  mean  molecular  weight.  The  high  resolution  of  the  laser 
device  offers  the  possibility'of  more  exact  temperature  profiles  (due 
to  small  angle  of  divergence). 

2. 1.4. 3  Measurement  of  Atmospheric  Constituents 

The  measurement  of  atmospheric  constituents  by  the  laser  seems  a 
good  possibility,  particularly  in  the  lower  atmosphere,  eg.  ozone: 
Unfortunately,  good  lasers  in  the  far  ultraviolet,  <  2500  A  are  not 
practical  but  in  the  near  ultraviolet,  ca.  3000  such  lasers  may  be 
developed.  .Also,  in  the  infrared  regions  many  irasers  will  be  available 
shortly'  Because  of  the  infrared  absorption  of  many  of  the  atmospheric 
gases  0^,  CO.^,  etc.,  these  irasers  should  be  effective  probes  The 


vertical  distribution  of  the  constituent  would  be  determined  from  the 
scattering  curve  which  would  not  behave  in  an  exponential  fashion  and 
hence  the ‘deviations  would  be  attributed  to  absorption.  Numerical 
analyses  have  been  made  of  this  possibility  for  the  carbon  arc  and  will 
not  be  repeated  here  since  the  computing  technique  is  similar. 

The  Russians  have  already  made  prototype  experiments  of  this  kind 
(without  laser)  using  UV  broadband  sources. 


2 , 2  The  Mass  Spectrometer 


Herzog  et  al .  of  GCA  have  constructed  and  flown  rugged ized  versions 
of  the  Redhead  mass  spectrometer.  In  essence,  the  instrument  permits 
the  measurement  of  the  profile  of  a  single  atmospheric  constituent.  The 
sensitivity' Of  the  instrument  is  such  as  to  permit  the  measurement  at 

atmospheric  pressure  of  an  impurity  of  one  part  in  10^^  or  at  much  lower 

4  3 

pressures,  a  density  of  10  particles/cm  . 


The  spectrometer  has  been  operated  extensively  in  the  laboratory 
and  its  performance  has  been  quite  satisfactory.  Figure  2-4  shows  a 
helium  peak  obtained  with  the  system.  Figure  2-5  shows  the  output 
currents  as  a  function  of  the  pressure  in  the  ion  source  The  slight 
deviation  of  the  linearity  is  caused  by  outgassing  and  the  remote  loca¬ 
tion  of  the  pressure  gauge,  and  will  be  eliminated  in  the  near  future 
Considerable  work  has  been  done  to  im.prove  the  electron  multiplier. 
Electrical  breakdown,  which  occurred  frequently  in  the  initial  experi¬ 
ments,  was  eliminated..  It  was  caused  by  marginal  dimensions  uf  the 
insulators  which  have  been  incorporated  in  the  first  flight  model. 
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OUTPUT  CURRENTS 

PRESSURE 

for 

helium-air  mixture 


PRESSURE  (tort) 


The  electronic  equipc.ent  associated  with  the  helium  mass  spectrometer 
may  be  conveniently  classified  into  four  main  categories:  the  electrometer 
amplifiers  and  the  three  sub-assemblies  that  comprise  the  main  power 
supplies,  a  regulated  26-volt  supply,  a  spectrometer  supply,  and  a  thermal 
ifagulator  system  to  eliminate  voltage  drift  due  to  temperature  changes. 

All  units  have  been  operated  satisfactorily. 

Minzner  has  pointed  out  that  on  the  basis  of  an  assumed  hydrostatic 
equilibrium,  the  temperature  profile  can  be  obtained  from  the  density 
profile  of  a  single  constituent.  Of  course,  the  sensitivity  of  such 
a  determination  would  depend  upon  the  scale  height  (and  therefore  the 
mass)  of  the  selected  specie. 


Total  densities  could  be  obtained  from  a  summation  of  several 
spectrometers  or  from  the  forenamed  laser  instrumentation. 

REFERENCES  ’ 

2-1  Schawlow,  A.  L,  and  Toxmes,  C.  H. ,  Phys .  Rev.  112.  1940  (1958). 

2-2  Javan,  A.,  Bennett,  W.  R. ,  Herriott,  D.  R. ,  Phys.  Rev.  Letters  6, 

106  (1961). 

2-3  Eltern-an,  L.,  "The  Measurement  of  Stratospheric  Density  Distribution 
with  the  Searchlight  Technique",  Geophysical  Research  Paper  No.  10, 
GRD  (Dec.  1951). 

2-4  Johnson,  E.  A.,  Meyer,  R.  C.,  Hopkins,  R.  E.  and  Mock,  W.  H,,  "The 
Measurement  of  Light  Scattered  by  the  Upper  Atmosphere  from  a 
Searchlight  Beam",  J.C.  S.A.  29,  512. 

2-5  Johnson,  E.  A.  and  Jcr.es,  RRDE,  Report  No,  367. 


2-6  Friedland,  S.  S.,  Katcenstein,  J.,  Sher.man,  J.,  Zatzick,  M.  R. , 

"Improved  Instrumenrat ion  for  Searchlight  Probing  of  the  Stratosph 
Physics  Dept.  Univ.  of  Conn,  Final  Research  Report  (.'■lay  195o). 


3i 


3. 


Solar  Flux,  Photo-Ionization  and  Absorption  Cross  Seccions 


The  photon  flux  incident  upon  a  region  of  the  upper  atmosphere  is 
determined  essentially  by  the  solar  flux  (outside  the  earth's  atmosphere), 
the  integrated  absorption  along  the  path  (as  determined  by  composition 
and  photo-absorption  cross  sections),  and  the  photo-ionization  cross  section 
in  the  region  of  interest.  These  parameters  then  determine  the  source 
function  for  ionization  (in  the  rf-breakdoun  problem)  which  is  external 
to  the  imposed  rf-field. 

The  gross  characteristics  of  line  emission  and  the  continuum  in  the 
solar  spectrum  have  been  delineated  in  the  review.  The  roles  of  solar 
UV  and  x-rays  have  been  motivated  on  a  semi-quantitative  basis.  In 
addition,  the  photo-ionization  and  absorption  cross  sections  have  been 
reviewed  and  listed  according  to  the  present  state-of-the-art. 

There  are  some  notable  deficiencies  in  our  information.  For  example, 
we  do  not  have  information  on; 

(1)  The  secular  and  spatial  variations  in  Che  solar  spectrum  outside 
of  the  visible. 

(2)  The  cornplete  absorption  spectra  of  Che  aCraospheric  species  below 

1  • 

luOu 

These  parameters  can  be  obtained.  Photon  counters  and  filters  for 
important  lines  such  as  monitoring  of  Lyman-'^  are  availab le ..  These  should 
be  used  on  satellites  to  provide  a  continuous  monitoring  of  solar  radiation. 
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The  absorption  work  can  be  done  In  the  laboratory  with  improved  light 
sources  But  the  question  is,  "Are  the  variations  in  these  parameters 
critical  to  the  rf-breakdoxm  problem?"  The  answer  at  present  is  that 
the  sophistication  of  the  theory  is  so  limited  that  these  parameters 
although  inportant  are  not  introduced  into  the  calculations.  Consequently 
it  would  appear  that  they  indeed  might  even  be  deferred  to  a  later  stage 
in  the  development  of  the  theory. 


4, 


Reaction  Rates 


Because  of  the  obvious  importance  of  the  field  of  reaction  rates  to 
the  behavior  of  the  atmosphere,  both  ambient  and  perturbed,  the  field  is 
being  very  actively  pursued  at  the  present  time.  The  existence  of  such 
a  broad  program  will  then  in^ose  certain  ground  rules  for  any  suggestions  . 
on  the  improvement  of  the  current  state-of-the-art.  In  short,  these  are; 

to  refrain  from  duplicating  existing  experiments  and  to  avoid  expensive 

*  • 

and  complicated  techniques  that  might  require  a  considerable  extension 
of  present  capabilities  and  therefore  involve  long  periods  of  time  before 
concrete  results  are  obtained.  Within  this  framework,  the  following  two¬ 
fold  program  involving  both  atmospheric  seeding  and  laboratory  instrumen¬ 
tation  is  suggested. 

4.1  Reaction  Rates  from  High  Altitude  Chemical*<Releases 

(4-1) 

Marmo  et  al  have  indicated  that  chemical  releases  in  the  upper 

atmosphere  under  varying  conditions  such  as  altitude,  degree  of  solar 
illumination  (using  Che  atmosphere  as  a  controllable  screen),  specie,  etc., 
can  be  used  effectively  to  get  rates  for  many  of  the  reactions  pertinent 
CO  the  dynamics  of  the  upper  atmosphere. 


In  a  ts'pical  ch 


(90-120  km)  alkali  m 
T  h  C  V  d  r  £  r  £  L  5  3  5  £  d  3  c 


mical  release  in  Che  altitude  regime  of  interest  here 
cals  such  as  sodium,  cesium,  or  lithium  may  be  used, 
high  temperature  which,  combined  with,  solar  radiation 
occ-sses,  results  in  the  creation  of  a  small  percentage 


of  free  electrons  The  chemicals  expand  initially  to  apprcximarp ? v  f-- 

1  ;  _  .  . 

atrrofpheric  density  and  then  continue  to  expand  radially  (but  less 

rapidly)  due  to  diffusion  processes.  During  this  latter  phase  of  the 

expansion  the  electron  population  is  depleted  by  the  various  attachment 

processes.  Radar  observations  of  the  expanding  electron  cloud  then 

measijre  the  rate  of  electron  depletion.  From  this,  from  the  altitude 

deperjdence  discussed  above,  from  the  fact  that  well-defined  chemical 

(or  Electron)  clouds  can  be  seeded  at  selected  altitudes,  and  from  the 

fact  ; that  a  wide  selection  of  suitable  chemicals  for  seeding  purposes 

are  at  the  disposal  of  the  investigator,  one  should  be  able  to  exliact 

reliable  values  for  chemical  reaction  rates  of  interest. 

In  general,  the  success  of  this  tool  depends  upon  several  factors; 

(1)  •  The  ability  to  produce  and/or  reproduce  a  particular  chemical 
release  (e.g.,  reaction  temperature,  completeness  jf  reaction,  etc,), 

(2)  The  ability  to  observe  significant  characteristics  of  the  chemi¬ 
cal  cloud  (electron  density  by  radar,  atoms  by  optical  emission,  etc.), 

(1)  The  ingenuity  of  the  researcher  in  the  fabrication  of  mcdsi? 
such  that  the  various  proces.ses,  say,  the  loss  processes  of  electrons 
(diffusion,  reccmbinat ion ,  attachment,  etc.)  can  be  readily  delineated 


These  three  vital  aspects  have  been  achieved  and  have 

of  evaluation  and  interpretation  of  both,  the  initial 
Cu-?) 

ny  Marmo  et  .ii'  and  the  much  more  inclusive  serie 


formed  the  basis 
series  of  releases 
3  of  releases  bv 
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Rosenberg  et  It  might  be  pointed  out  thet  the  conclusicn  of  t^is 

series  of  tests  does  not  mean  that  there  has  been  full  utilization  of  the 
electron  cloud  as  a  tool  for  investigating  the  parameters  of  the  upper 
atmosphere.  There  are  other  uses  (of  the  contaminant  release)  which  have 
provided  the  motivation  and  predominated  in  the  design  of  the  tests. 

(»  .  * 

The  phenomenology  of  the  release  depends  upon  such  parameters  as 
chemical  yield  of  the  contaminant,  thermal  ionization  efficiency, 
ambipolar  and  neutral  diffusion,  photo-ionization  cross  sections,  chemical 
consumpclon,  etc.  Depending  upon  the  pertinence  of  the  interaction  to 
a  particular  release,  the  investigator  must  be  able  to  pull  out  piecemeal 
each  of  these  parameters,  It  would  be  remiss  to  enter  into  a  detailed 
description  of  each  and  the  appropriate  solution  of  each  type  of  release 
as  a  function  of  altitude.  Actually,  the  dynamics  of  the  release  is  given 
by  the  continuity  equation,  and  the  solutions  of  a  score  of  these 
characteristic  equations  were  given  in  Chapter  7  of  Atmospheric  Processes. 
However,  it  is  advantageous  to  demonstrate  that  each  atmospheric  para¬ 
meter  can  be  distilled  out  if  sufficient  ingenuity  is  applied.  Consider, 
for  example,  the  parameter  of  diffusion.  Figure  (^-1)  shows  the  appli¬ 
cation  cf  a  particuiar  mathematical  model  to  a  series  of  releases  which 
varied  in  altitude.  Note  that  the  selected  mathematical  model  of  diffu¬ 
sion  control  is  completely  successful  in  giving  the  time  history  of  the 
release  for  release  altitudes  of  IIG  to  133  km  Ac  9i  and  100  km,  ocher 
interactions  such  as  chemical  consumption  and  attachment  are  important 
enough  to  materially  effect  the  electron  density.  But,  and  this  is  the 
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essence  of  the  entire  natter,  the  releases  under  a  diffusion-controlled 
configuration  have  made  it  possible  to  separate  out  the  diffusion  contri¬ 
bution  in  the  complex  release.  In  a  similar  manner,  each  and  every  atmos¬ 
pheric  parameter  including  rates  of  pertinent  chemical  reactions  can 
be  defined  or  evaluated. 


4,1.1  Radiative  Attachment  to  Atomic  and  Molecular  Oxygen 

An  examination  of  a  typical  ^reaction  is  offered* as  an  example  of  the 
dynamics  of  contaminant  releases  in  the  upper  atmosphere,  the  choice  being 
attachment  to  oxygen.  The  dominant  processes  of  attachment  of  electrons 
in  the  pertinent  region  of  the  upper  atmosphere  are:’ 


two-body  radiative  attachment 
0  r  e  *♦  0  +  h  v 

O2  +  e  O2  hv 

three-  body  attachment 


Bioiidi  well  S'?  Rloch  and  Rradhury  at  an  earlier  dace)  has  measured 

in  Che  Laboratory  the  chree-faody  race  of  attachment  for  molecular  oxygen, 
the  resulting  race  being  of  the  order  of  3  x  IG  cm^  sec  Similiar 
laboratory  measurements  on  the  races  of  radiative  attachment  in  molecular 
oxygen  have  not  been  achieved  in  the  laboratory.  One  reason,  of  course, 
is  the  relative  speed  of  the  competitive  process  of  three-body  attachment. 
In  order  to  circumvent  this  diificulty,  it  is  necessary  Co  go  to  much 


» 

'•  'J' 


smaller  pressures,  but  then  the  radiative  flux  becomes  too  small  to 
measure  and/or  wall  effects  become  Important.  As  for  atomic  oxygen, 
measurement  on  such  a  specie  in  the  laboratory  is  not  realistic  since 
the  lifetime  is  too  small. 

How  radiative  attachment  rates  of  both  atomic  and  molecular  oxygen 
can  be  accomplished  by  cogent  experiments  in  the  upper  atmosphere  is 
sketched  as  follows. 

]^en  a  point-release  electron  cloud  is  formed  the  chemicals  expand 

almost  Immediately  to  a  particle  density  approximately  that  of  the  ^bient 

after  which  they  diffuse  outward  much  more  slowly.  Man  release  experiments 

have  indicated  that  the  radial  particle  distribution  in  the  electron  cloud 

is  approximately  Gaussian  and  remains  so  for  a  relatively  long  period  of 

time.  The  important  parameter  for  this  distribution  is  the  Gaussian 

•• 

half-width  r^.  The  second  important  parameter  is  the  diffusion  velocity 
which  is  implicitly  contained  in  the  diffusion  coefficient  D=  Both  r 

-  ■  o 

and  D  have  been  determined  rather  accurately  for  many  different  altitudes 
and  for  many  different  types  of  chemical  clouds.  Above  an  altitude  of 
about  100-110  km  (see  Figure  4-1)  the  chemical  cloud  expansion  is  clearly 
diffusion  controlled  and  chemical  consumption  and  turbulent  processes  are 
unimportant.  However,  below  these  altitudes  the  latter  two  processes  do 
become  important,  and  though  complicate  the  analysis,  Che  contributions 
of  chase,  processes  are  separable. 

The  essential  characteristics  of  releases  in  the  exponentially  varying  - 
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VALUES  FOR  SEVERAL  RELEASES 


(- 


DIFFUSION  MODEL  OOO  EXPERIMENTAL 


Figure  4-1 


»  Q 


atnospheres  are: 

(1)  In  the  altitude  regine  between  about  90  km  and  120  km  the 
oxygen  content  in  the  atmosphere  suffers  a  marked  change.  At  90  km  the 
concentration  of  molecular  oxygen  is  about  100  times  that  of  atomic 
oxygen.  As  the  altitude  increases  above  90  km  the  molecular  oxygen 
photo-dissociates  so  that  at  about  105  km  the  concentrations  of 
molecular  and  atomic  oxygen  are  approximately  the  same,  and  at  120  km 
the  concentration  of  atomic  oxygen  is  about  10  times  that  of  molecular 
oxygen.  These  facts  are  clearly  depicted  in  Figure  4-2  where  the 
atmospheric  concentrations  of  atomic  and  molecular  oxygen  are  shoiim  as 
functions  of. altitude.  The 'importance  of  this  profile  is  that  it  enables 
the  investigator  to  sort  out  three  general  altitude  regimes;  chemical 
Interactions  Involving  molecular  oxygen  predominate  say  between  90  and_ 
perhaps  100  km,  those  involving  atomic  oxygen  predominate  above  115-120  km, 
and  interactions  involving  both  atomic  and  molecular  oxygen  are  jointly 
important  for  altitudes  like  100-110  km. 

(2)  Atmospheric  temperature  and  pressure  both  vary  considerably 
between  90  and  120  ktn.  The  pressure  dependence  should  provide  the 
means  of  isolating  three-body  reactions  (these  are  more  important  at 
the  higher  pressures,  hence  at  the  lower  altitudes)  from  the  two-body 
radiative  reactions  which  should  be  mors  import. snt  at  the  higher 
altitudes.  To  what  extent  chemical  release  measurements  can  be  refined 
so  as  CO  yieid  temperature  data  remains  to  be  dsterm.ined. 


lao 
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(3)  Both  accniic  and  molecular  oxygen  have  strong  positive  electron 
affinities  whereas  molecular  and  atonic  nitrogen  have  negative  affinities 

(4)  Atmospheric  constituents  other  than  oxygen  and  nitrogen  are 

present  in  only  very  small  amounts.  • 

(5)  In  the  altitude  regime  between  about  90  and  120  km  chemical 
seeding  techniques  are  particularly  amenable  to  chemical  reaction  rate 
analysis.  This  is  because  chemical  consumption  is  minimal  and  because 
diffusion  rates  are  relatively  slow  thus  allowing  for  long  observation 
times. 

■  Item  (1)  should  provide  the  means  of  isolating  the  rates  of  radiative 
attachment  to  atomic  and  molecular  oxygen.  Item  (2)  should  enable  one 
to  distinguish  two*body  radiative  attachment  from  three-body  attachment. 
Item  (3)  can  distinguish  reactions  involving  oxygen  from  those  involving 
nitrogen.  Item  (4)  extends  this  distinction  to  all  constituents  other  - 
than  oxygen.  And  Item  (5)  extends  the  analysis  to  constituents  which  do 
not  naturally  occur  in  the  upper  atmosphere  yet  retains  a  correct 
environment  free  of  wall  effects. 

It  is  possible  to  state  the  case  in  a  sophisticated  manner;  for 
measurements  at  a  sufficiently  high  altitude  so  that  threa-body  attach¬ 
ment  is  not  important  but  sufficiently  low  so  that  the  presence  of 
molecular  oxygen  is  still  important  the  attachment  probability  is 

p(h)  =  kjol  +  [O2]  (4-1) 
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where  the  brackets  [  ]  indicate  concentrations.  Then  from  several 
measurements  of  p  at  different  altitudes,  and  from  a  knowledge  of  [o] 
and  [O^l,  we  can  in  principle  determine  kj^  and  k^. 


At  the  lower  altitudes  tdiere  three-body  attachment  predominates 
slmiliar  expressions  for  rate  constants  of  three-body  attachment  to 
atomic  and  molecular  oxygen  could  be  obtained. 


Collision  cross  sections  may  be  simply  obtained  from  ihe  rate 


constants  in  case  the  are  temperature  and  pressure  independent. 
,  ifie  following  formula  is  then  generally  valid: 


k.  =  a.  V.  e 
1  i  ir 


D./kT  A  ■  N 


(4-2) 


where  is  the  interaction  energy  and  is  the  relative  velocity. 


4.1.2  Associative  Detachment 


A  second  example  of  the  scheme  is  the  reaction  of  associative 
detachment 


0  s 
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Dalgarno  states  that  for  this  reaction  a  rate  as  high  as  lO"*’*^  cn^  sec"^ 
cannot  be  e.scluded.  Yet  a  siniple  examination  of  the  data  on  nigh  altitudi: 


chemical  releases  indicates  that  a  maximum  rate  for  this  reaction  is 


10  cm  sec 


6',. 

For  an  opClmuin  effort  in  rates  by  high  altitude  releases,  new 
eKperiments  would  be  conducted,  and  perhaps  even  new  tools  such  as  the 
Puerto  Rico  dish  (employing  incoherent  back  scatter  of  electrons) 
applied.  However,  in  the  modest  effort  proposed  here,  the  costs  of  con¬ 
ducting  supplementary  tests  are  circumvented. 

4 . 2  Reaction  Rates  from  Laboratory  Program 

Marmo  has  suggested  a  group  of  experiments  to  measure  some  of  the 
present  unknown  rates  employing  such  selective  techniques  as  photo- 
dissociation.  For  example,  consider  the  much  discussed  reaction  of  atom- 

ion  exchange,  which  was  originally  assigned  a  rate  by  Bates  as  high  as 

-9  3  -1  -12  3  -1 

10  cm  sec  (later  estimated  as  10  cm  s  ) 

0“^  +  NO"^  +  0 

It  has  been  customary  in  measuring  rates  of  such  reactions  to  create 
the  ion  specie  by  electron  impact.  Unfortunately,  such  a  means  is  not 
particularly  selective,  and  a  number  of  ionic  species  and  excited  states 
are  formed  in  addition  to  the  one  desired.  This  complicates  the  measure¬ 
ment  considerably.  Marmo’s  idea  is  to  employ  specific  photo-disscciation 
and  photo-ionization  to  create  a  particular  specie.  In  the  above  measure 
-ent,  commencing  with  a  mixture  of  N2  and  O2,  we  have: 

O2  ^  hv  (1470  or  8.45  ev)  -  0  (^D)  *  0  (^P) 

0  ('’0)  -  hv  (11.6  ev)  -»  0  m  e 

0*-  -*•  +  N 


The  ionization  and  dissociation  energies  of  the  gases  involved  are; 


0(^P) 

13.6 

0(^D) 

11.6 

^2 

12.2 

5.12 

”2 

15.6 

9.76 

NO 

9.-^ 

6.5 


The  energies  of  the  selected  photons  are  therefore  such  that  formation 
of  NO  proceeds  by  a  unique  process.  The  ratio  [NO^]  [O]  as  measured 
in  a  mass  spectrometer  can  then  be  related  to  the  rate. 

In  summary,  then,  a  modest  program  combining  simplifed  laboratory 
measurements  with  the  detailed  study  of  release  of  contaminants  in  the 
upper  atmosphere  gives  definitive  results  which  can  improve  the  existent 
state-of-the-art  for  reaction  rates. 
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5.  Transition  Probabilities 

In  the  sequence  of  events  leading  to  brea!«down  by  the  coupling  of 
rf  energy  to  a  slightly  ionized  gas,  inelastic  collisions  between  elec¬ 
trons  and  neutrals  result  in  the  excitation  of  electronic,  vibrational, 
and  rotational  levels.  The  lifetimes  of  these  levels  (essentially  the 
inverse  of  the  transition  probability)  then  determine  the  energy  loss 
due  to  radiative  transitions.  Whether  or  not  the  radiated  photon  es¬ 
capes  from  the  medium  depends  upon  the  population  of  the  lower  (energy) 
state  and  the  photo-transition  probability.  The  relative  (oscillator) 
strengths  for  the  various  means  of  excitation  are  given  roughly  by  the 

characteristic  figures  for  electronic  excitation,  0.1  to  0,5  ev  for 

-2  -4 

vibrational  levels,  and  10  to  10  ev  for  rotational  levels.  These 

figures  are  to  be  compared  to  the  average  energy  gain  for  electrons 

between  collisions;  essentially  all  the  energy  is  randomized  in  a 

single  collision,  the  factor  being  (1-m/M).  In  air  at  an  equivalent 

7  2 

altitude  of  30  km,  the  breakdown  field  is  of  the  order  of  10  watts/m  . 

2 

Dividing  by  the  flow  velocity  yields  the  energy  density  (€  E  )  which 

3 

is  of  the  order  of  O.Oi  joules/m  ,  or  E  ~  400  volts/cm.  The  mean  free 

-4 

path  of  electrons  at  this  altitude  is  about  5  x  10  cm,  so  that  the 
energy  gain  per  collision  is  of  the  order  of  10  "  volts.  Since  this 
is  of  the  same  order  of  this  vibrational  energy  increm.enc,  the  impor- 
t.ance  of  the  cross  sections  for  these  inelastic  collisions  is  clearly 
dem.onscrated .  It  appears  that  electronic  excitation  in  oxygen  occurs 
at  2.7  ev,  but  there  is  also  evidence  of  vibrational  excitation  at  a 
few  tenths  of  an  ev .  The  manner  in  which  these  various  factors  enter 
into  the  Boltznwinn  equation  is  given  in  Section  7.2. 
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The  subject  of  transition  probabilities  of  all  atomic  species  which 
are  relavent  to  the  ambient  atmosphere  (an  atmospheres  contaminated  by 
particular  chemical  seedings)  has  been  adequately  treated  by  A.  Naqvl  in 
Atmospheric  Processes.  To  give  an  estimate  of  the  accuracies  "conferred" 
upon  the  values  for  atomic  transitions,  it  is  considered  that  the  theore¬ 
tical  values  are  good  to  +  10%  while  the  experimental  uncertainty  is  + 
30%.  In  spite  of  the  Inability  of 'the  experimentalists  to  substantiate 
the  theoretical  results,  the  latter  are  considered  to  be  on  firm  grounds 
established  by  the  quantum  theory. 

The  question  now  is  -  what  is  to  be  done  to  advance  the  state-of- 
the-art  so  as  to  provide  information  On  transition  rates  which  may  be 
essential  to  a  complete  understanding  of  the  mechanics  of  breakdown. 

There  are  three  steps  which  suggest  themselves  and  they  are  given  in 

the  order  of  increasing  comple.xity  as  well  as  degree  of  importance  to 

\ 

the  breakdo'.jn  problem.  These  are:  First,  to  extend  the  calculations 
on  atomic  species  to  other  hea’.-y  atoms.  These  atoms  can  result  from 
nuclear  detonations  but  the  extension  may  be  more  of  esthetic  than  prac¬ 
tical  value.  The  second  and  more  important  step  is  to  calculate  the 
transition  probabilities  of  the  molecular  species  (e.g.,  0^,  N^,  NO,  etc.) 
Tlie  former  can  be  accomplished  by  a  relatively  simple  extension  of  the 
screening  theories  of  Layzer  (Section  5.1).  However,  the  problems  intro¬ 
duced  by  going  to  molecular  species  (Section  5.2)  are  inherently  more 
complicated.  The  third  step  is  to  obtain  the  cross  sections  for  excita¬ 
tion  of  tliese  various  states  by  electron  impact.  Unfortunately,  there 
is  no  unambiguous  procedure  of  ennloying  the  results  of  the  more  tract¬ 
able  photo-cxcitation  to  the  process  of  excitation  by  electron  impact. 


The  probability  an^litude  for  a  radiative  transition  A  -*  B  in  a 
nany-electron  atom  is  proportional,  in  a  first  approximation,  to  the 
matrix  element  ,  where  D  is  the  total  electric  moment 

of  the  atom  and  u  is  the  polarization  vector  of  the  absorbed  or  emit¬ 
ted  photon.  If  onjs  makes  the  further  approximation  of  assigning  the 
states  A  and  B  to  definite  configurations,  the  matrix  element  reduces 
to  a  product  of  two  factors,  one  involving  spin  and  angular  coordinates 
only,  the  other  involving  radial  coordinates  only.  Since  all  matrix 
elements  connecting  a  given  pair  of  configurations  have  the  same  radial 
factor,  the  spin-angular  factor,  which  can  be  evaluated  by  means  of 
standard  methods  in  the  theory  of  atomic  spectra,  determines  the  rela¬ 
tive  strengths  of  Zeeman  components  in  a  line,  in  lines  in  a  multiplet, 
and  of  multiplets  in  a  transition  array.  Thus  the  determination  of 
absolute  transition  probabilities  hinges  on  the  radial  factor. 


One  can  allow  for  configuration  interaction  by  assigning  the 
states  A  and  B  to  mixtures  of  pure  configurations.  Once  the  components 
in  such  a  mixture  have  been  chosen,  their  relative  amplitudes  can  be 
evaluated  with  little  difficulty.  The  transition  matrix  element  ■ 

<^'u  -  may  therefore  be  regarded  as  a  known  linear  combination 

of  matrix  elements  connecting  pure  configurations,  and  the  problem  re¬ 
duces  formally  to  one  in  which  configuration  interaction  no  longer 
figures  explicitly,  though  in  fact,  as  we  shall  see,  its  most  important 

from  Davi.j  t_iy;:cr  (iXV;) 
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effect  still  remains  to  be  taken  into  account. 


The  radial  factor  is  proportional  to  the  transition  integral 


Kni.  n'J')  =  r  P„.^.dr  , 


<3-l) 


where  F  ,  and  P  .  are  the  radial  wave  functions  of  the  active  electron 
ni  n  i 

in  States  A  and  B  respectively.  The  factor  of  proportionality  is  a 
product  of  integrals  of  the  type  * 


I  (nf)  =  P^,  P®,  dr, 

o  nf  nf 

o 


(5-2) 


one  for  each  inactive  electron,  and  is  normally  very  close  to  unity. 


Although  the  radial  transition  integral  depends  explicitly  on  the 
radial  functions  of  a  single  electron,  each  of  these  depends  implicitly 
on  the  entire  set  of  radial  functions  for  the  state  in  question.  Thus, 
if  configuration  interaction  is  neglected,  the  best  radial  functions  for 
a  given  atomic  state  satisfy  a  set  of  simultaneous  integrodifferential 
equations  (the  Fock  equations),  a  typical  member  of  which  has  the  form 


r . 


Z  (Z-i-l) 


z  -  s  (r) 

IlA/ 


2  2 
dr 


P  ,  =  E  ,  P  ,  i-  X  (5-3) 
n*  n^  n^  nZ 


The  screening  function  S  ,,  and  the  exchange  function  X  ,  are  actually 
functionals  of  all  the  radial  functions  in  the  set,  and  their  forms  de¬ 
pend  on  the  quantum  numbers  S  and  L  as  well  as  on  the  configuration 
quantum  numbers. 


Allou-ir.g  for  configuration  interaction  by  assigning  the  states  A 


and  3  to  specific  mixtures  of  configurations,  one  arrives  by  way  of  the 


variation  principle  at  a  set  of  "generalised  Foe k  equations".  A  typical 
member  of  this  set  has  exactly  the  same  form  as  a  typical  Fock  equation, 
but  the  exchange  term  in  the  generalized  equation  Includes  contribu¬ 
tions  arising  from  configuration  interaction  as  veil  as  from  exchange. 
This  is  the  effect  of  configuration  interaction. 

Numerical  solution  of  the  Fock  equations,  even  in  their  original 

form,  presents  too  many  practical  difficulties  to  serve  as  an  adequate 

/  . 

basis  for  extensive  calculations  of  f-values.  Besides,  such  a  procedure 
vould  tell  very  little  about  the  accuracy  of  the  results.  One  needs  to 
simplify  the  equations,  not  only  in  order  to  be  able  to  solve  them,  but 
also  in  order  to  gain  insight  into  how  various  aspects  of  their  structure 
affect  the  accuracy  of  transition  integrals  evaluated  from  their  solu¬ 
tions  . 

Owing  to  the  form  of  the  integrand,  the  value  of  a  typical  transi¬ 
tion  integral  depends  strongly  on  the  behavior  of  the  radial  functions 
at  large  distances  from  the  nucleus.  This  behavior  is  in  turn  largely 
deLeniiined  by  the  energy  parameters  relating  to  the  active  electron  in 
the  initial  and  final  states.  Finally,  the  energy  paratneters  are  related 
to  ionization  energies  whose  values  can  he  found  by  experiment.  All 
these  relations  are  of  fundamental  importance  to  the  problem. 

5.1.1  Asymptotic  Form  of  the  Radial  Functions 

As  r  increases,  S  ,  —  N-l,  where  N  is  the  number  of  electrons  in  the 
n<J 

system,  and  X  ,  0,  so  chat  the  Fock  ecuacion  for  P  .  reduces  to  the 

ni  *  ni 

radial  wave  equation  for  an  electron  in  a  Coulomb  field. 
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1  d"  iliii)  Z  - 

2  2^  2 

dr  2r 


I 

S  ->-1  I 

J 


ni  n2  n2 


(5-4) 


It  then  follows  that 


(n*  ^  zy. 


w  M 

,  n*,i-l/2  V  n  / 


(5-5) 


(Z  =  Z  -  N  *  1) 


where  t^\l2  IJhittaker's  hypergeometric  function,  and  n*,  the  effec¬ 


tive  principal  quantwtz  number,  is  defined  by 


(5-6) 


The  normalization  factor  A  has  been  chosen  so  as  to  have  the  value  I 
when  n*  is  an  integer.  Apart  from  this  factor,  the  asymptotic  form  of 
thus  depends- entirely  on  the  value  of  (or  n*). 

The  parameter  has  a  very  simple  interpretation.  Let  denote 
a  given  electron  configuration  minus  an  nZ  electron,  so  that  (U)  =  (U  ) 
(nf) .  IJhen  the  systen  as  a  whole  is  in  a  definite  energy  eigenstate  F, 
the  subsystem  whose  configuration  is  will  not  normally  be  in  a  definite 
one  of  its  energy  eigenstates  but  in  a  mixture  of  them.  Let  E(U  jp)  de¬ 
note  the  expectation  value  of  the  energy  of  the  subsystem  when  the  system 
as  a  whole  is  in  the  state  F.  If  E(U!r)  denotes  the  absolute  energy  of 

I  , 

the  State  F,  then  £  _  -  EfUjF)  -  E(U  iF).  Thus  E  ,  is  the  expectation 
value  of  the  energy  of  the  system  relative  to  the  "core"  (defined  as  the 
system  minus  an  nj  electron).  Tiiis  result  remains  valid  if  the  state  F 
is  assigned  to  a  mixture  of  configurations  instead  of  just  one. 


The  energy  of  the  core  is  not  the  same  as  the  energy  of  the  ion 
formed  by  removing  an  nl  electron  from  the  atom.  The  former  is  always 
numerically  greater  (smaller  in  magnitude)  than  the  latter,  because  the 
nZ  electron  tends  to  shield  its  companions  from  the  field  of  the  nucleus. 

It  is  convenient  to  define  the  energy  of  a  system  relative  to  its 
parent  ion,  E*(U^|r),  by  analogy  with  the  definition  of  E(U^|r).  Thus 
E*(U*jr  )  is  minus  the  expectation  value  of  the  minimum  energy  needed  to 
ionize  the  system.  Though  unconventionally  defined,  this  quantity  does 
admit  in  principle  of  direct  experimental  determination.  In  practice, 
only  the  energies  of  the  ionic  states  are  known,  not  their  probabilities. 
However,  the  latter  can  be  calculated  by  means  of  standard  methods  in 
the  theory  of  complex  spectra.  We  may  therefore  look  upon  the  energy 
of  a  system  relative  to  its  parent  ion  as  an  observable  quantity.  De¬ 
noting  this  quantity  by  we  have 

j 

w^^  =  E(Ulr)  -  E*(U^ Id,  (5-7) 

so  that 

w  ,  -  E  „  =  E(u"lr)  -  E*(u'^!r)  >  o.  (5-8) 

n  A/  rii/ 

Tai s  di  fference  represents  the  effect  of  screening  by  the  ni’  elec¬ 
tron  on  the  energy  of  the  core  and  need  not  be  small  compared  with  -W  ,, 

n^ 

Some  recent  quantitative  studies  show  that  in  neutral  atoms  fW  ,  -  E  .) 

n2  nf’ 

is  apnreciable  compared  with  -W  .whenever  the  core  contains  one  or 

no 

more  electrons  that  belong  to  the  sa.me  shell  as  the  screening  electron, 
and  increases  i.’ith  the  number  of  such  electrons.  For  the  ground  state 


of  ;:0,  with  its  eight  L  electrons,  the  contribotion  of  2p  electron  to 
the  core  energy  turns  out  to  be  several  tit^s  as  great  as  the  ionization 
energy  itself.  These  studies  also  show  that  the  relative  iaportance  of 
the  screening  contribution  diminishes  with  increasing  degree  of  ioniza¬ 
tion  along  an  Isoelec tronic  sequence. 

When  the  principal  quantum  number  of  the  active  electron  exceeds 
the  principal  quantum  numbers  of  all  the  core  electrons,  its  contribu¬ 
tion  to  the  core  energy  is  usually  small,  though  not  always  negligible. 

5.1.2  The  Coulomb  Approximation 

A  transition,  and  the  states  it  connects,  is  "qussi-hyarcgeriic"  if,  in 
borh  states,  the  principal  quantum  number  of  the  active  electron  exceeds 
the^ principal  quantum  numbers  of  all  the  core  electrons.  Except  in 
alkali-like  atoms,  transitions  to  or  from  the  ground  states  of  many- 
electron  atoms  rarely  qualify  for  this  label.  On  the  other  hand,  quasi- 
hydrogenic  transitions  between  excited  states  are  very  common. 

^fast,  chough  probably  not  all,  quasi -hydrcgenic  transitions  satisfy 
the  following  two  conditions; 

I,  Tlie  active  electron  shields  the  cere  electrons  so  ineffec¬ 
tively  that  W  ,  =5=  E  ,  and  W  ,  ,,  =5  £  . 

n  c  n  i  n  '  ^  n '  ^  * 

2:  The  value  of  the  transition  integral  depends  almost  entirely 

on  the  contributions  from  a  range  (r  .oo)  within  which  the  radial  func- 

C  ' 

Cions  cay  be  placed  by  their  asymptotic  forms. 


If  these  conditions  are  fulfilled  and  if  the  energies  W  ,  and  W  ,  , 

tiv  n  *  -w  * 

are  knoun  from  experiments,  the  transition  integral  can  be  evaluated, 
apart  from  the  normalization  factors  A  and  A',  since  it  can  be  approxi¬ 
mated  by  an  integral  over  the  range  (r^,  a>),  within  which  the  integrand 
is  knoim  to  a  good  apprcrtimation.  Extensive  tables  of  such  Integrals 

have  been  published  by  Bates  and  Damgaardf^  these  incorporate  earlier 

(5-2) 

calculations  of  the  satae  kind  by  %lleraas.  The  principles  under¬ 

lying  the  approximation,  however,  have  been  exploited  much  earlier  by  a 
number  of  authors,  though  in  a  less  explicit  and  less  systematic  fashion. 


Let  r  .  denote  the  average  nuclear  distance  of  an  electron  whose 
ni  or 

radial  function  is  P  ,(r  .=  P  ,rP  .dr),  and  let  p  denote  the  geometric 
ni  ni  ni  nZ 

o 

mean  of  r  ,  and  r  j .  The  pure  number  1  -  Q  with  Q  =  I(n^,  n'i’)/p, 
measures  the  extent  of  cancellation  between  positive  and  negative  contri¬ 
butions  to  the  transition  integral  I.  By  Schwarz's  lemma,  Q  <  1,  with 
equality  only  when  nJ  *  n‘2*.  The  smaller  Q  is,  the  more  stringent  be¬ 
comes  the  second  of  the  two  conditions  given  above. 


This  condition  may  he  written  in  the  form 


P  ,  -  ?  .  ..  dr  «  Q. 
n .  c  n  i  • 


(5-9) 


If  Q  1.  this  inequality  is  weaker  chan  the  conditions 


?  T ..  dr  «  1 


(5-10) 


o 


o 


for  the  nornalization  factors  A  and  A’  to  be  close  to  1,  so  that  in  some 
circumstances  one  may  be  able  to  Improve  the  approximation  by  using  more 
accurate  values  for  A  and  A'  than  A  =  A"  =  1.  Seaton'  has  shown  how 

to  derive  such  values  from  the  energy  spectra  in  which  and  ^^e 

imbedded . 

■  .  -  '  -*1 

Even  when  the  two  conditions  that  are  necessary  for  the  validity  of 
the  Coulomb  approximation  are  satisfied,  the  quantum  defects  n-n*  and 
n’-n‘*  may  be  substantial.  Only  for  very  highly  excited  states  are  the 
quantum  defects  negligible.  These  states  are  hydrogenic,  not  merely  quasi 
hydrogenlc .  For  them,  the  Coulomb  approximation  is  trivial.  Thus  the  do¬ 
main  within  which  the  approximation  can  usefully  be  applied  seems  to  have 
rather  clear-cut  boundaries. 

5,1.3  Beyond  the  Coulomb  Approximation 

Before  turning  to  specific  proposals  for  going  beyond  the  Coulomb 
approximation,  let  us  consider  in  a  general  way  what  this  would  involve. 

Transitions  in  which  the  active  electron  belongs  to  the  same  shell 
as  one  or  more  of  the  inactive  electrons  do  not  normally  meet  either 
of  th.e  two  conditions  previously  mentioned.  If  the  active  electron  be¬ 
longs  to  the  same  subshe  1 1  as  one  or  m.ore  of  the  inactive  electrons,  it 
is  almost  certain  chat  neither  is  fulfilled.  This  means  that  we  can  no 

longer  rely  on  exaeriment  to  furnish  E  ,,  and  therefore  we  can  at  most 

n  w 

hope  to  develop  methods  that  Will  apply  to  transitions  for  wh.ich  the 
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energy  parareters  can  be  predicted  with  sufficient  accuracy."^  Moreover, 
the  effects  of  the  nonconstant  part  of  the  screening  function  and  of 
the  exchange  function  on  the  form  of  the  racial  functions  can  no 
longer  be  ignored ■ 

Because  one  can  always  evaluate  energy  parameters  more  simply  and 

with  greater  precision  than  wave  functions,  it  is  useful  to  regard  the 

evaluation  of  E  '  as  the  first  step  toward  a  determination  of  P 
n^  ■  n£ 

Given  E  ,  together  with  suitable  approximations  to  S  .  and  X  ,  one 
n£  nf  nz 

could  integrate  the  equation  for  inward  from  infinity.  The  solu¬ 
tion  would  deteriorate,  of  course,  as  it  approached  the  origin,  and  ulti¬ 
mately  diverge;  but  such  solutions  could  still  yield  accurate  estimates 
for  transition  integrals. 

*• 

5.1.4  Application  of  the  Screening  Theory 

As  is  well  known,  the  predictions  of  the  standard  theory  of  atotaic 
spectra  concerning  absolute  term  energies  are  subject  to  large  errors 
of  a  systetnatic,  nature ,  arising  from  the  neglect  of  configuration  inter¬ 
action.  Recently,  an  alternative  theoretical  scheme  has  been  proposed 

rs ./. ) 

for  describing  complex  spectra  (Layzer ''  ),  whose  principal  features 

are  as  follows: 

1.  By  allowing  the  nuclear  charge  7.  to  enter  explicitly  into 

^’’Sufficient  accuracy"  means  in  practice  that  the  effective  principal 
quantum  numbers  n'-  and  n’*  must  be  known  to  within  small  fractions  of 
unity.  The  smaller  the  value  of  Q,  the  more  accurately  the  energy 
par-ameters  need  to  be  known. 
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all  the  calculations,  one  is  able  to  classify  the  effects  of  configuration 
interaction  in  a  physically  significant  manner  and  to  work  consistently 
to  a  given  order  of  approximation. 

2.  Certain  interactions  between  configurations  are  taken  into 
account  in  the  first  approximation,  namely,  those  connecting  states  be¬ 
longing  to  the  same  "complex",  a  complex  being  defined  as  the  aggregate 
of  all  states  belonging  to  a  given  set  of  principal  quantum  numbers  and 
a  given  parity. 

Within  the  framework  of  the  approximation  defined  by  the  require¬ 
ments  that  Z  enter  explicitly  into  all  the  calculations  and  that  all 
interactions  connecting  terms  in  the  same  complex  be  taken  into  account, 
the  simplest  theory  is  that  based  on  screened  Ig^drogenic  radial  functions. 
The  best  screening  parameters  are  those  defined  by  the  variation  principle 
In  spite  of  its  simplicity,  the  screening  theory  accounts  quantitatively 
for  the  main  regularities  characterizing  term  spectra  in  isoelec tronic 
sequences.  Its  predictions  are  also  free  from  the  systematic  errors 
associated  with  the  usual  theory. 

Tlie  predicted  tern  energies  have  the  form 

W  =  W  Z~  -  W,  Z  ^  W  ^  0(z"S.  (5-11) 

Z  n  O 

apart  from  relativistic  corrections.  Tae  first  two  coefficients  are 
given  e.xactly  by  the  theory,  the  rentlining  ones  only  approxirnately.  Pne 
absolute  accuracy  of  these  predictions,  so  far  as  it  has  been  tested, 
appears  to  be  quite  good  in  s«derately  and  highly  ionized  systems. 
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In  neutral  atons,  the  terms  of  order  1/Z  are  often  conspicuous,  and  the 
agreement  between  the  screening  theory  and  experiment  is  less  good. 

Let  us  now  return  to  the  problem  of  determining  suitable  radial 
functions  for  evaluating  transition  integrals.  To  begin  with,  we  are  now 
in  a  position  to  fill  a  gap  that  we  left  earlier  in  our  discussion  Of  the 
generalized  Fock  equations.  Ve  assumed  that  the  states  A  and  B  could 
each  be  assigned  to  definite  mixtures  of  pure  configurations.  It  is  now 
clear  that  we  must  allow  for  all  interactions  within  a  given  complex. 

The  mixture  coefficients,  or  relative  amplitudes  of  the  configurations 
in  a  given  mixture,  are  then  independent  of  Z  and  their  values  are 
given  by  the  screening  theory.  We  do  not  in  this  way  allow  for  all  the 
effects  of  configuration  interaction,  but  we  allow  for  the  most  impor¬ 
tant  ones  and  we  ensure  the  correctness  of  the  first  two  terms  in  the 
expansions  Of  the 

In  a  first  approximation  One  could  use  values  for  E  ,  and  E  , ,, 

Hi  n  xi 

calculated  according  to  the  screening  theory.  As  for  the  troublesome 
terms  S  ,  and  X  ,,  these  are  of  relative  order  1/Z  compared  with  the 
other  terms  in  the  equation  for  ?  .,  Probably  it  would  be  a  sufficiently 

n^r 

good  approximation  to  use  screened  hydrogenic  functions  to  evaluate  them. 

Ihe  equation  for  ?  ,  would  tl;ca  reduce  to  either  a  single,  independent, 
n 

linear  integro-dif fercr.tial  equation  or  to  an  inhomogeneous  differential 

equation,  denending  on  how  one  dealt  with  the  occurrences  of  P  ,  in  the 

nc 

functional  .  In  either  case  a  numerical  solution  would  present  no 
difficulty  The  nuclear  charge  Z  would  appear  e.xplicitly  in  the  equation. 
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and  could  probably  also  be  retained  in  the  solution.  Thus  a  single 
numerical  solution  would  furnish  values  of  transition  integrals  for  homol¬ 
ogous  transitions  in  an  entire  isoelectronic  sequence. 

5.1.5  Transition  Integrals  in  Highly  Ionized  Systems 

In  an  approximation  that  becomes  increasingly  accurate  as  Z  increases 
along  an  isoelectronic  sequence,  one  can  evaluate  transition  Integrals  by 
u^ing  screened  hydrogenic  radial  functions,  the  screening  parameters  being 
given  by  the  screening  theory.  Vaxsavsl^^^  has  carried  out  a  detailed 
comparison  between  the  results  obtained  with  this  method  and  those  ob¬ 
tained  by  means  of  more  elaborate  and  (presumably)  more  accurate  methods. 
For  certain  classes  of  transitions  (e.g.,  transitions  in  which  the  prin¬ 
cipal  quantum  number  undergoes  no  change)  he  found  a  remarkably  close 
correspondence  between  the  two  sets  of  predictions,  extending  even  to 
transitions  in  neutral  atoms. 

Since  the  comparison  data  were  largely  of  unknown  accuracy,  the 
meaning  of  this  result  is  not  entirely  clear.  However,  to  some  e.xtent 
one  can  estimate  the  accuracy  of  transition  integrals  calculated  on  the 
screening  approximation  independently.  A  given  integral  can  easily  be 
evaluated  as  an  explicit  function  of  Z  and  of  the  screening  parameters 
for  the  initial  and  final  states.  Tnis  integral  is  accurate  to  terms  of 
relative  order  \.ii-  Nov;,  small  changes  of  the  screening  parameters  will 
also  produce  changes  of  relative  order  l/Z  in  the  value  of  the  integral  . 
Moreover,  one  can  estimate  the  inherent  uncertainty  in  the  values  of  the 
screening  parameters,  Pnus  one  can  estimate  one  of  the  contributions  to 


the  error;  and  it  is  not  unreasonable  to  suppose  that  the  renaiaing  con¬ 
tributions  will  usually  be  of  the  same  order  of  magnitude.  In  this  way 
we  arrive  at  a  simple  and  useful  criterion:  If  the  value  of  a  transition 
integral  is  insensitive  to  small  changes  in  the  screening  parameters  it 
is  fairly  reliable. 

5.2  Transition  Probabilities  of  Molecular  Species 

In  Atmospheric  Processes,  the  Morse  potential  functions  for  0^,  N^, 
and  NO  as  calculated  by  Forrest  Gilmore  (Rand  Corporation)  were  presented. 
It  was  indicated  (e.g.,  by  inclusion  of  the  repulsive  state  in  )w0  suggested 
by  Marmo)  that  the  listing  of  the  Morse  curves  were  incomplete.  Recently, 
Miescher^^  has  discovered  several  more  of  these  repulsive  states. 

These  should  be  explicitly  added  to  the  overall  picture  since  their  inclu¬ 
sion  is  essential  to  the  proper  interpretation  of  experiments  on  lifetime 
of  metastable  state,  etc.  Gilmore  also  listed  approximate  values  of 
these  lifetimes  (inverse  of  the  transition  probability)  of  some  of  the 
metastable  states  involved.  Again,  these  listings  are  both  incomplete 
and  insufficiently  well  determined.  In  the  interactions  of  waves  and 
plasmas,  in  addition  to  the  electronic  transition  probabilities,  one  re¬ 
quires  the  vibrational  and  rotational  transition  probabilities 

Tne  theoretical  vibration  eigenvalues  E  can  be  derived  from  the 

V 

series : 

(v  r  1/2)  -  h  =  X  (v  -  1/2)- 
e  e  e 

♦  h  c  .  y  fv  -  I /2 .  (5-12) 


E 

V 


ii  C 


c 


which  essentially  deterisines  the  shape  cf  the  potential  well. 

Usually  and  w  are  experimentally  determined.  The  vi¬ 
brational  wave  functions  are,  of  course,  interrelated  to  the  vibrational 
eigenvalues  and  the  assumed  shape  of  the  potential  function,  that  of  the 
Morse  function  being  # 


V(r)  =  1 


-a  (r  -  r  ) 
e  e 


(5-13) 


The  Morse  potential  leads  to  an  exact  solution  of  the  wave  equation 
H  terms  of  Laguerre  polynomials,  but  the  Franck-Condon  terms 

converge  so  slowly  that  the  calculation  of  the  vibra- 

V  V  j. 

2  -  2 
tional  transition  probabilities  R^,^„  =  dr  becomes  very 

cumbersome.  In  the  above,  R  is  the  electronic  transition  moment. 

e 

Among  the  various  suggestions,  Bouigue^^  Fraser  et  al^^ 

Bates and  the  suggestion  of  Wu  concerning  the  possible 

application  of  the  WKB  (phase  integral)r.ethod  may  provide  a  convenient 
solution.  Hie  WKB  method  has  the  advantage  of  being  valid  from  any 
arbitrary  potential  function  which  is  not  as  radically  asymetric  or  has 


large  gradients. 


It  also  autoroaticallv  takes  into  account  the  shift 


of  the  intemuclear  distance  (f^)-  Furry has  indicated  chat  in 
the  zeroth  vibrational  level,  the  deviation  of  the  WKB  amplitudes  are 
of  the  order  of  2\  For  higher  vibrational  levels,  the  agreement  should 
gradually  improve 


The  vibrational  quantum  numbers  of  the  oscillator  having  Che  poten¬ 
tial  function  Vfv)  arc-  given  by  the  phase  integral 


where  the  integrand  is  the  classical  momentum  of  the  oscillator,  and 
Xj  and  are  the  classical  turning  points  (the  momentum  vanishes). 

Forms  of  V(x)  have  been  given  by  Hulbert  and  Hirschfelder 

The  WKB  method  then  gives  the  approximate  solution  to  the  vibra¬ 
tional  wave  function  as  exponential  in  the  two  regions  outside  of  the 
turning  points  and  a  cosine  function  inside  these  limits,  From  the 
wave  functions,  the  vibration  transition  probabilities  as  given  approxi¬ 
mately  by 

®J’v"  “  ‘**^1  (5-15) 

can  be  obtained  by  computer  methods. 

One  may  conclude  that  it  will  be  no  simple  task  to  get  all  the 
photo  and  (electron)  impact  transition  probabilities  required  in  the 
Boltzmann  equation  for  the  interaction  of  an  rf  source  x^ith  a  slightly 
ionized  gas.  Yet,  this  .must  be  accomplished  before  this  mathematical 
discipline,  x.’hich  can  in  pri.nciple  describe  the  co.T.plicated  state, 
is  fruitful . 

5.3  Excitation  by  Electron  I.mnacf-' 

Tile  quantum  mechanical  solution  of  electron  impact  x.;ich  an  ato.mic 
specie  has  received  considerable  attencioa  -  at  least  for  the  case  of  the 
hydrogen  atom.  Tee  hamilconian  for  the  tv;o  electrons  involved  is  given  by 


=  (v  +  1/2)  jr 


(5-14) 


=  .  .  t!  f',72  _  el 

2n  v  1  2)  ^2 


(5-16) 


To  obtain  solutions  to  the  wave  equation  (H^  =  e'{'  ),  it  is  customary  to 


utilize  an  approximation  of  the  form 


^2>  =  ^/)  ’^n  ^^1>  <^2> 


(5-17) 


wherein  the  degeneracy  between  impact'-  and  bound  electron  is  conveniently 


ignored .  The 


(‘  •  f ) 


takes  into  account  both  discrete  and  continuous 


states  respectively.  Substituting  the  approximate  form  of  in  the  wave 
equation,  multiplying  by  the  complex  conjugate  and  integrating  over 
space,  one  obtains  the  infinite  set  of  simultaneous  equations 


where 


(v-^  -t-  k  =  (e  +  f  )  U  F 

\2  n/n  ^  J  J  nmm 

„  2me^  /I  1  \  , 

U  =  'if*  ( -  -  —  14  dr, 

nm  .1.2  .  n  Vi,  ^2'  “  ^ 


(5-18) 


2 

k  ^  (£  -  €  ) 

n  h  n 


Asymacotic  solutions  of  the  subsidirry  functions  (F  )  can  then  be  ob- 

n 

tained  frcn:  the  usual  criteria  that  the  wave  function  is  well  behaved, 

-and  that  the  interaction  vanishes  for  large  distance  (r,);  these  have  the 


f  ’  t  )  -  e 


i^o-r 


(5-19) 


'Subscript  2  denotes  tb.e  irnact  electron. 


For  the  solid  angle  increaent  d'J)  (0,~),  the  differencial  scattering  cross 

‘2 

section  for  excitation  of  state  i  is  I  =  If  (5,5)1  . 

n  '  n 

For  fast  impact  electrons  (i.e.,  uhen  the  energy  of  the  impact  elec* 

tron  is  much  greater  than  the  binding  energy  of  the  bound  electron)  the 

perturbation  is  small,  so  that  the  Born  approximation  is  valid.  However, 

for  the  case  of  rf -breakdown,  fast  electrons  are  not  involved,  and  the 

Born  approximation  is  inappropriate.  A  weaker  restriction  than  that  of 

the  small  perturbation  (the  two-state  approximation  with  close  coupling) 

is  to  assijote  that  except  for  those  terms  associated  with  the  initial 

atomic  state,  the  diagonal  terms  of  the  matrix  U  predominate  over  the 

nn 

non -diagonal  ones.  Physically,  this  implies  that  the  distortion  of  the 
initial  and  final  "wave"  of  the  impact  electron  by  the  mean  static  fields 
of  the  atom  in  both  its  initial  and  final  state  is  included,  but  that 
intermediate  states  of  the  system  are  ignored.  Under  these  conditions, 
the  wave  function  and  the  set  of  differential  equations  become: 

-  "o'  ■  “oo)  ^0-0 

-  Q  ’  “„o  fo 

^  -  Vq  (-^)  Fq  (r^)  -  (r^)  (r^)  (:;-20) 

CoiT'-'T.enc ing  from  the  equations  can  be  solved  by  iteration,  che  dis- 
Cor*ad  \;ave  being  given  by  the  first  tern  in  the  series. 

The  degeneracy  between  the  impact  and  bound  eiec crons  cam  be  taken 
into  account  in  the  t'..'o-state  approximation  by  imcluding  terms  which  are 

qT  - 


« 


either  sy::s=ecric  or  antisycnetric  (denoted  by  superscript  «ith  respect 
to  spatial  interchange.  The  ’/  function  and  the  wave  equation  then  yield: 


± 


^  -  =  (r^)  (r^)  =  (r^)  (r^)  +  (r^)  F^"  (r2) 


(5-21) 


+  Icq  -  Uqq^  Fq  (r^)  ±  Kqq  (r^  F^  (r^  d  r* 


U  F  "  ^ 
nO  n 


•^On  ^  ^2 


(5-22) 


^  n  nnj  n  2  nn  2  2a  2  2 


=  «n0^0-  ^  (^2  'P  ^0‘  (^2>  ^  ^2 


(5-23) 


where  K„„,  K„  and  K  are  given  by 
00  On  nn 


V  =  (r.)  V*  (r  )  jk  ~  -  ^  —  +• 

inn  n  1  n  2  n  OJ 


^nrn  (5-24) 


In  keeping  with  the  two-state  concept,  and  K  represent  additional 


distortions  of  the  initial  and  final  electron  ’..-aves,  while  K.  is  a  coupling 

Un 


contribution  which  induces  transition  between  the  tw-o  states.  Tne  differ¬ 
ential  cross  section  now  takas  the  form 


I  (},  i)  da  - 
n 


Pf 


(^)r 


f  (^): 


d.j 


(5-25) 


In  this  exchange  and  distorted  wave  approximation,  wheie  and  are 


small,  the  iteration  proceeds  rapidly  yielding,  say  for  the  process  of 
elastic  scattering 


-  “oo)  ±  "^0  '2>  ^0'  <'2>  ■'  'I  ■  “ 


th 


and  for  inelastic  scattering  (for  n  state) 


*  k  ^  -  U  ^  F  -  +  K 
\  n  nny  n  “  i 


nn  V 


=  “nO  ^O"  ^nO  ^2^  ^1  '^2 


(5-27) 


Bates and  Seaton^^  have  shovm  that  in  the  case  of  colli¬ 


sions  of  slow  electrons  v;ith  atomic  oxygen,  additional  terms  must  be  - 
inserted  to  allow  for  close  coupling  to  Intermediate  states  of  Che  same 
configuration.  If,  further,  the  small  energy  difference  between  these 
states  is  ignored,  the  problem  remains  a  tractable  one. 


The  restriction  appropriate  to  the  several  approximations  can  be 
listed  as: 

Distorted  wave  -  =  0;Uq^  is  small 


Fotential  distorted 


exchansc 


So  "  Sn  ■-  Sn’  Sn 


Exchange  and  potential 


distorted  exchange  -  are  small 


Close  Coupling  uichouc 
e::chan^e 


■  ^00  ’■  \n  “  Sn  “  ° 


Close  coupling  with 
exchange 
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From  an  examination  of  the  results  of  the  various  disciplines  for 
the  relatively  simple  case  of  the  hydrogen  atom,  it  appears  that  the 
distorted  two-state  approximation  is  sufficiently  accurate,  provided 
that  the  exchange  interaction  is  included.  For  transitions  involving 
several  states  of  nearly  the  same  energy,  coupling  between  states  must 
be  considered. 


Insofar  as  the  rf -breakdown  is  concerned,  the  interest  is  quite 
naturally  in  atomic  and  molecular  species  such  as  N^,  O^,  NO,  0 
rather  than  atomic  hydrogen.  For  these  cases,  it  is  more  realistic  to 
resort  to  measurement  to  obtain  the  necessary  cross  sections  for  ex¬ 
citation. 


In  "Atmospheric  Processes",  cross  sections  of  only  a  few  of  these 
pertinent  electron-neutral  excitations  (or  de-excitations)  were  listed. 
Tliese  Included: 

From  Seaton^^  ^^^(for  atomic  oxygen) 

^P  -  *S  ; 

for  atomic  nitrogen 


■?  4 

D  -  S 


From  Perceval  (for  atomic  ozrygen) 

3  3 

P  •  3a  P 


P  -  3p  P 


From  Stewart 

0  -  0,  0  -  1,  and  0-2  bands 

^2  (X-  Lp  (bV) 

This  brief  summation  gives  some  idea  of  the  dearth  of  data  in  the 
areas  of  electronic,  rotational,  and  vibrational  excitation  of  perti¬ 
nent  species  (especially  N2,  O2  and  the  oxides)  under- electron  impact. 
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6 ,  Radiation  Losses  In  RF  Excited  Plasmas 

The  rate  of  energy  loss  by  radiation  is  essential  to  the  determin¬ 
ation  of  breakdown  threshholds.  In  particular,  there  are  at  least  two 
types  of  radiative  processes  to  be  considered.  The  first  and  dominant 
process  includes  all  microscopic  interactions  such  as  atomic  and  molecular 
radiative  transitions  and  perhaps  to  some  extent  chemiluminescence.  The 
theoretical  determination  of  the  transition  probabilities  was  discussed 
in  Section  5.  Experimental  determinations  of  radiative  transitions  in 
an  excited  gas  whose  ionization  varies  over  a  wide  range  are  extremely 
difficult.  Monitoring  of  the  decay  of  microwave  discharges  is  being  done 
but  data  are  meager.  The  basic  difficulty  is  that  at  the  low  densities 
of  the  upper  atmosphere  where  radiation  can  become  more  ia^sortant  than 
the  collisional  processes,  the  observations  must  be  made  over  long 
optical  paths.  In  the  laboratory,  the  3-body  (2-body  plus  wall)  effects 
mask  the  slower  two-body  radiative  processes. 

Upper  atmosphere  experiments  on  radiative  loss  processes  are  in 
general  very  expensive  when  the  investigation  provides  the  energy  (rf) 
of  excitation,  However,  several  of  these  experiments  have  already  been 
funded;  chat  they  are  not  equipped  to  give  the  maximum  information  on 
optical  effects  is  unfortunate.  The  reason  is  simply  that  the  dominent 
motivation  of  these  exper Lii'.ents  is  enhanced  electron  density  and  its 
decay  modes.  Of  course,  natural  air-glow  and  aurora  can  give  some 
information,  but  here  excitation  is  chiefly  from  electron  impact,  and  the 
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distribution  of  excited  states  may  vary  widely  from  that  of  a  stratum 
irradiated  by  an  rf-beam  because  of  the  wide  variations  in  the  electron 
energy  distributions. 

•  • 

The  three  high  power  density  experiments  currently  in  (near)  opera¬ 
tion  are: 

(1)  The  Bailey  experiments  in  Australia  which  beams  .5  megawatts 

of  power  at  the  gyro-frequency  (about  1.5  mc/s)  vertically  incident  at 
the  ionosphere.  The  basic  idea  is  simply  to  raise  the  electron  temperature 
of  a  thin  stratum  absorption  at  the  gyro- frequency  being  particularly 
effective.  If  three-body  attachment  (e  +  0^  "^02^  important  electron 

loss  process,’  the  electron  density  will  increase  (Bailey  calculates  a 
factor  of  10)  since  the  attachment  cross  section  decreases  with  energy 

in  this  range.  The  principle  "observation"  will  be  the  electron  density 
as  m.easured  by  a  C-4  ionosonde. 

(2)  At  Jicamarca,  Peru,  Ken  Bowles  is  conducting  high  power  density 

I 

irradiation  of  th.e  upper  atmosphere  ^t  a  frequency  of  50  mc/s  and  a  power 
of  six  megawatts.  Evidently,  NBS  is  aware  of  the  possible  importance  of 
the  enhanced  "air-glow".  However,  and  this  is  an  important  point,  the 
tests  are  at  present  unrelated  to  the  breakdowTi  problem.  For  this  rea¬ 
son,  there  will  be  no  attempt  to  optimize  the  equipment  for  the  purposes 
suggested  herein.  Furthermore,  equipment  for  the  several  experiments  will, 
if  no  action  is  taken,  be  completely  different  for  each  installation,  and 
no  steps  will  be  taken  to  utilize  the  results  of  all  the  experiments  in 
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furthering  the  furdaaiental  understanding  of  the  breakdown  phenomenon. 


(3)  In  Puerto  Rico,  a  1000  ft  "dish”  is  being  constructed  for  studies 
on  incoherent  back  scatter  from  electrons  (not  fluctuations  in  refractive 
index)  at  a  frequency  of  400  mc/s.  Again  high  power  densities  are  avail¬ 
able;  in  particular,  the  focusing  of  the  beam  is  much  better  than  that 
obtainable  in  the  ether  experiments.  However,  the  amount  of  absorption 
may  be  significantly  decreased  since  the  non-deviative  absorption  coeffi¬ 
cient  (k)  is  given  5y 


2n 

me  4  (i  ui 


foJ, 


(6-1) 


Only  in  the  lower  atmosphere  are  electron  collision  frequencies  high 

enough  (v  >  u)  to  make  k  basically  independent  of  u.  Otherwise  for 

2 

uj  >  V,  k  varies  as  1/cj  ,  and  the  absorption  of  the  rf  will  be  very  small 
-  -  perhaps  well  beyond  detection  of  radiative  transitions.  The  Puerto 
Rico  experiment  is  therefore  not  an  optimum  one  for  our  purposes.  Yet, 
if  no  action  be  taken  here,  no  attempt  will  be  made  to  employ  photometric 
devices  to  monitor  the  experiments  in  the  eventuality  that  optical  infor¬ 
mation  is  available. 


In  addition 
composition  and  s 
is  described  by  t' 
inhomogeneities  1; 
which  may  affect 


Che  microscopic  character  determined  essentially  by 
e  occ'jpat  Lons ,  a  plasma  under  an  external  Ell  field 
macroscopic  field  vectors.  Under  these  imposed  fields, 
he  plasma  can  give  rise  to  energy  loss  by  radiation 
electron  density  at  and  near  breakdown.  As  evidenced 


1 
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£a  the  review,  some  theoretical  work  has  been  done  on  the  problem,  but 
this  work  has  not  been  conclusive  in  establishing  the  importance  of  this 
mechanism  to  the  breakdown  phenomenon.  Therefore,  it  is  suggested  that 
a  theoretical  study  of  the  problem  be  made,  commencing  with  a  review 
and  evaluation  of  past  work.  It  seems  likely  that  such  action  can  place 
the  mechanism  in  proper  perspefctive. 


»  O 
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7.  Electron  Energy  Distribution  Functions  and  the  Boltzmann  Equation 

The  time  seqtieiice  of  events  leading  to  r£  breakdown  can  only  be 
fully  described  by  the  Boltzmann  Equation.  Unfortunately  for  real 
cases,  a  total  solution  of  the  problem  which  must  include  all  the 
microscopic  and  macroscopic  interactions  is  beyond  oiir  present  con^uta- 
tional  capabilities  and  perhaps  even  our  knowledge  of  the  specific 
microscopic  interactions  involved.  It  is  quite  natural  then  that  the 
*.  total  problem  is  greatly  simplified  by  theoretical  workers  who  make  an 

effort  to  phrase  the  problem  in  such  a  manner  that  the  problem  is  tract¬ 
able.  The  question  raised  by  these  efforts  (a  subject  often  inadequately 
stated  in  the  literature)  is:  what  is  the  motivation  and  range  of  validity 
of  the  various  approximations?  In  the  remainder  of  Cection  7,  R.  Papa  has 
reviewed  the  various  theoretical  treatments,  and  delineation  and/or 
evaluated  their  relation  to  the  problem  of  rf  breakdown.  This  material 
supplements  Chapter  XIII  (Plasma  Dynamics)  in  Atmospheric  Processes. 

7.1  The  Mechanics  of  RF-3reakdo^TO 

If  an  electromagnetic  field  is  impressed  upon  a  partially  ionized 
gas,  the  free  charges  will  gain  energy  from  the  EM  field  by  suffering 
collisions  with  neutral  species  and  having  their  ordered  oscillatory 
motion  changed  to  random  thermal  motion.  In  the  lower  atmosphere,  the 
free  charges  will  exist  initially  net  as  electron-ion  pairs;  but  as  small 
ions,  Aitken  muclei,  and  molecular  ions  (no  thermal  electrors  exist  very 
long  since  their  lifetime  at  STP  is  10  sec  due  to  three  body  attachment 
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rate  of  5  x  10  cm^  sec  ^).  These  free  aolecular  ions  will  eventually 

give  rise  to  electrons  through  Inelastic  collisions  with  neutral  species 

(collisions  with  neutrals  in  xnetastable  states  may  be  an  effective 

electron  production  mechanism) .  In  the  upper  atmosphere,  the  daytime 

3  3 

ambient  electron  density  varies  from  10  particles  per  cm  in  the  D 

5  3 

region  to  10  particles  per  cm  in  the  F  region. 

As  the  electron  temperature  increases  due  to  randomization  of  Che 
electron  oscillatory  motion  through  elastic  collisions  with  neutrals, 
there  will  be  various  electron  energy  loss  processes  such  as: 

(1)  Collisions  which  can  excite  electronic  levels  of  atoms  and/or 
molecules 

(2)  Collisions  which  excite  molecular  vibrational  levels  (sub¬ 
excitation  electrons) 

(3)  Collisions  which  excite  molecular  rotational  levels  (subexcita¬ 
tion  electrons) 

(4)  Ionization  of  atoms  and  dissociative  or  non-dissociative 
ionization  of  molecules. 


Processes  in  itec;  (4)  actually  increase  the  aunber  of  electrons  present  and  a 
knowledge  of  the  energy  distribution  function  of  the  secondaries  produced  by  the 
primaries  is  necessary  for  a  detailed  calculation  of  the  breakdown  condition. 

The  electron  concentration  may  also  increase  through  collisional  and  associative 
detachs^nt.  In  Che  upper  atmosphere,  photodetachment  is  an  additional  electron 
production  mechanism.  Electron  loss  processes  include  dissociative,  radiative  and 
three-body  recombination,  and  radiative  and  three-body  attachment.  A  microscopic 
description  of  the  interaction  of  an  EM  field  with  a  gas  (in  terms  of  Che  electron 
energy  distribution  function)  involves  a  knowledge  of  the  cross  sections  for  the 
following  processes:  elastic  electron-neutral  collisions,  ionization  and  excitation 
collisions, radiative  and  collisional  deactivation  of  excited  states,  radiative,  three- 
body  and  dissociative  attachment,  electron-ion,  ion-ion  and  three-body  recombination, 
charge  transfer,  ion-atom  exchange,  and  collisional  and  associative  detachment. 

Even  if  a  thorough  knowledge  of  the  various  rates  and  cross  sections  as  a  function 
of  energy  permitted  an  elimination  of  some  of  the  relatively  unimportant  processes, 
the  formulation  of  the  problem  in  terms  of  a  Boltzmann  equation  for  the  electron 
energy  distribution  function  is  extremely  complicated  and  includes  the  following 
effects: 

(1)  In  cha  anbienC  atmosphere  there  are  various  macroscopic  processes  such  as 
turbulent  diffusion.,  winds,  shears,  and  fluid  transport  forces  on  a  region  where  the 
E3-I  radiation  is  b-eing  absorbed. 

(2)  The  neutral  composition  may  change  (by  such  processes  as  N2  +  e  -*  N2'*'  +  2e 
and  then  e  +  N.,'  —  N  +  N)  even  before  breakdown  is  reached. 


Most  of  the  expericuital  work  on  RF-breakJown  has  consisted  of  studies  with 
simple  geometrical  arrangem.ents  of  EM  fields  in  a  microwave  cavity.  In  which  the 
breakdown  condition  be  expressed  as  a  straightforward  boundary  value  problem. 
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The  gases  e’r.plcyed  in  the  laboratory  experiments  are  usually  either  hydrogen  or 
pure  air,  and  a  radioactive  source  is  used  as  the  external  ionizing  agency  in  place 
of  the  solar  flux.  In  most  of  the  laboratory  experiments,  electron  loss  by  re¬ 
combination  may  be  neglected  compared  with  loss  by  attachment  and  diffusion.  Brown 
and  MacDonald^^  have  formulated  the  CW  RF-breakdown  condition  for  a  particular 
gas  for  a  specific  EM  mode  in  a  microwave  cavity  in  terms  of  (1)  the  characteristic 
value  of  the  time  Independent  continuity  equation  for  the  electron  density 

7^  Ng  +  (U/D)  Nq  *  0 

where  V  =  u.  -  u 

i  a  .  . 

=  Frequency  of  ionization  due  to  electron  impact. 

=  Frequency  of  attachment  for  electrons  to  neutrals. 

D  =  Electron  diffusion  coefficient, 
and  (2)  the  high  frequency  ionization  coefficient 

C  =  P/DE^ 

where  .E  =  electric  field  anqilitude. 

The  CW  breakdown  condition  (neglecting  recombination  compared  with  attachment  and 
diffusion)  for  the  case  of  infinite  parallel  plates  with  a  uniform  electric  field 
may  be  expressed  as 

v/D  =  l/\- 


or 


where  A  =  L/.t  =  diffusion  length 
L  =  pla  te  separation. 

This  breakdown  condition  is  valid  within  the  following  three  limits: 


(1)  The  mean  free  path  is  less  than  any  dimension  of  the  cavity, 

(2)  The  frequency  is  sufficiently  high  so  that  the  electrons  do  net  lose 
appreciable  energy  between  cycles,  and 

(3)  The  average  motion  of  the  electrons  resulting  from  the  action,  of  the  field 
and  of  the  collisions  is  sufficiently  small  so  that  the  field  does  not  sweep  the 
electrons  out  of  part  of  the  cavity  in  each  half  cycle. 


One  of  the  most  general  conditions  for  SF  breakdown  may  be  described  by  the 

condition  dn  /dt  —  0  where  n  is  the  electron  density  (n  satisfied  the  continuity 
e  e  ■  e 

(7-2) 

equation).  Hartman  has  investigated  the  validity  of  such  an  approach.  In  his 

study,  a  theoretical  expression  for  the  electron  density  was  obtained  as  a  function 
of  the  RF  field  amplitude.  However,  Hartman  made  the  following  simplifying  assumptions: 

(1)  Ionization  occurs  as  a  result  of  single  impacts  between  gas  atoms 
and  sufficiently  fast  electrons. 

(2)  The  gas  has  an  infinite  volume,  and 

(3)  Negative  ions  are  not  formed,  so  that  the  only  mechanism  for  electron 
removal  is  recombination  with  positive  ions. 

His  calculations  show  that,  at  a  certain  critical  field  strength  (the  breakdown  field), 
the  density  of  electrons  rises  sharply.  His  simple  model  leads  to  a  breakdown 

phenom.enon  af  field  strengths  not  far  from  the  observed  values.  Th^  gases  treated 
are  helium  and  neon  at  frequencies  of  3000  megacycles.  In  Che  amb-ienC  atmosphere, 

the  condition  dn  /dt  0  will  also  describe  the  Cr<’  RF  breakdown  crirerion.  However, 

e  ■■■■,, 
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for  pulsed  RF,  the  brcakdo'..-n  criterion  nay  not  be  expressed  as  the  point  at  which 

the  rate  of  gain  of  electrons  (irpact  ionization  and  photo^ionization)  equals  the 

rate  of  loss  of  electrons.  For  a  pulsed  RF  beam,  once  breakdown  has  occurred  the 

electron  concentration  will  increase  at  the  rate  e'^^,  where  v  =  v.  -  V  -  D/A^ 

i  A 

and  =  frequency  of  ionization  by  electron  Impact 

V  =  frequency  of  attachment  of  electrons  to  neutrals 

D  =  electron  diffusion  coefficient 
A  =  diffusion  length  ’  ’ 

and  photoionization  and  recombination  have  been  neglected.  This  buildup  will  occur 
for  any  pulse  whose  amplitude  is  equal  to  or  greater  than  that  of  the  CH  breakdown 
field.  The  RF  pulse  can  penetrate  the  region  until  the  electron  density  has  reached 
the  plasma  resonant  density 

2  2 

n  =  m  <u  An  e 
•  p  e 

where  oi  =  frequency  of  EM  wave. 

Xt  may  be  concluded  from  the  previous  discussion  that  a  knowledge  of  the  time-dependent 
electron  energy  distribution  function  is  sufficient  to  determine  the  breakdown  con¬ 
dition  for  CW  and  pulsed  RF  power.  The  OJ  breakdo^m  condition  may  be  found  by  deter¬ 
mining  the  electric  field  arnplituca  for  which  the  electron  density 

r  3 

n^Ct)  =  J  f^Cv.t)  d  V 

shov.'o  a  sharp  increase  in  magnituca.  In  the  case  of  RF  pulses,  a  large  part  of  the 
R?  power  may  be  transmitted  throwsh  the  hot  spot  in  the  ambient  atmosphere  even 
though  breakdov.-n  has  occurred.  Hare  the  RF  pulse  will  be  reflected  when  the 
electron  density  reaches  the  plasma  resonant  density: 


82 


n 


f  (v.t)  d  V 
e 


n 

P 


n  oo“/4:te 
e 


The  determination  of  the  electron  energy  (or  velocity)  distribution  function  of 
course  involves  a  solution  of  the  Boltzmann  equation. 
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7.2  The  Boltziaan  Equation 


In  the  case  of  a  partially  ionized  plasma,  an  adequate  description  of  the 
energy  distribution  function  of  each  species  involves  the  solution  of  a  series 
of  coupled  Boltzmann  equations: 


e  ■* 


^+vVf  -  —  [E  +  vxB]-Vf  =  B  + 
ot  r  e  m  V  e  eo 

e 


I  'e  !«>.  ^  I 

A  J  If 


“k 


ion'£ 


n-1 


=  B 


e  Ion, 


+  y  B,  +  Y  E 

ion  io.L  Ion  n 

k=l  ^  ^  J  ■ 


df 


N-1 


+  V'V  f  =  Vb  ,  +B  +Vb 

at  r  nj  L  Hj  nj 


.th 


where  (f.  ),  is  the  distribution  function  for  the  I  ion  constituent,  f  is  the 

Ion  I 


.th 


distribution  function  for  the  j  neutral  constituent,  the  subscript  e  refers  to 


electrons,  subscript  n  to  neutrals,  and  B  is  the  contribution  to  the  collision 

mn 


th 


integral  for  the  distribution  function  of  the  m  type  of  particles  In  collision  with 
th 


an  n  type.  The  electron,  ion,  and  neutral  equations  are  coupled  through  the 
collision  Integral  terms,  3^^*  The  macroscopic  fields  which  appear  in  these 
equations  are  self  consistent.  In  that  the  electrons  are  not  conceived  of  as  inter¬ 
acting  directly  with  one  another.  Rather,  they  produce  an  electromagnetic  field 
(the  field  vectors  E  and  B  must  satisfyMar-rwell's  equations)  which  is  e:tpressad  as 
a  function  of  position  and  time.  The  electromagnetic  field  then  acts  on  each 
electron  through  its  position  and  tima  coordinates.  In  order  that  such  an  approach 


In  order  tha 


ba  valid,  it  is  necessary  that  the  mactoscopic  fields  E  and  B  reraain  constant  during 
a  collision  (ion-electron,  electron-neutral,  ion-=eutral).  It  is  further  assumed 
that  three-body  forces  are  negligible  compared  vita  two-body  forces. 


Most  of  the  papers  on  the  kinetic  theory  of  microwave  discharges  in  gases 

(Holstein^^”^\  Margenau^^  Allis  and  Brown^^  Reder  and  Brown^^  neglect 
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ion-electron  and  electron-electron  collisions.  However,  Ginzburg  and  Gurevich  , 
whose  work  is  butlined  in  sections  7.4  and  7.5  have  discussed  the  conditions  under 
which  it  is  permissible  to  neglect  electron-electroa  and  electron-ion  conditions 
in  a  partially  ionized  plasma.  Processes  such  as  photo-ionization  and  photo-detachment 
have  never  been  considered  in  the  literature  as  terms  in  the  Boltzmann  equation. 
Electron  elastic  collisions  with  neutrals  are  the  most  easily  handled  process. 
Expressions  for  the  collision  integral  for  clectroa-neutral  collisions  may  be  ob¬ 
tained  in  terms  of  the  neutral  density,  temperature,  elastic  collision  cross  section, 
and  zero  order  electron  energy  distribution  function.  The  derivation  of  these  expres¬ 
sions  usually  involves  the  assumptions  that  the  electron  suffers  only  a  small  change 
in  speed  on  collision  with  a  neutral  (basically  tha  factor  is  1-m/m) .  More  sophistLcated 

(7_8) 

treatments  (such  as  those  of  Allis'  ,  which  is  presented  in  Section  7.4)  take  into 
account  the  recoil  of  the  heavy  molecule  under  elastic  collisions  with  electrons. 

Almost  all  treatments  on  the  interactio.n  of  micro-vsve  energy  with  a  gas  consider 


only  the  changes  in  the  electron  distribution  function,  where  the  ion  and  neutral 


distribution  functions  are  assumed  to  remain  Maxwellian. 


Inelastic  collisions  between  electrons  and  neutrals  are  accompanied  by  the 
excitation  of  rotational,  vibrational  or  electronic  levels  and  also  by  ionization.  In 
addition,  second-order  impacts  are  possible,  in  which  the  energy  of  the  excited  state 
of  Che  molecule  is  transferred  to  the  incoming  electron.  Ho'..’ever,  an  exact  calculation 
of  all  these  inelastic  processes  would  be  very  ccmplicated.  The  appropriate  cross 
sections  are  known  only  in  a  few 


cases. 


7.3  Spherical  Harmonic  Expansion s 


The  preseace  of  the  collision  integral  terns  in  the  Boltzmann  equation 

(which  cakes  the  value  of  f  at  v  depend  on  the  number  of  electrons  at  very  different 

velocities)  characterizes  the  equation  as  a  non-linear  partial  differeatio-integral 

equation,  for  sviiich  no  general  method  of  solution  is  known.  It  is  customary  to 

employ  a  perturbation  technique,  by  expanding  the  function  f,  to  obtain  solutions 

(7-9) 

to  this  equation.  Enskog's  method  consists  in  expanding  the  distribution 

function  in  powers  of  a,  where  the  system  is  considered  to  be  close  to  thermal 
equilibrium  and  departures  from  equilibrium  are  caused  by  some  agent  a: 


12  2 

f  =  f  Q  +  af  +  aV  + . 


Unfortunately,  this  technique  does  not  converge  well  for  charged  particles  in  an 
“electric  field. 

For  the  case  of  electrons  (but  not  ions)  in  an  electric  field,  a  spherical 
harmonic  expansion  of  the  distribution  function  does  converge  more  rapidly: 

f(r,v)  =  y  f^  ir.v)  P.  (cos  9) 


.0  ,  .1 

=  f  T  f  cos  - 


where  the  electric  field  is  taken  to  lie  along  the  z-axis. 


Two  scalar  quantities  that  may  be  obtained  from  the  zero  order  distribution 
function  are:- 


T  .  .  ....  (7-8) 

Parts  ot  osettan  3.  are  ta.-ten  tro.a  Allis 


O 


the  electron  concentration: 


Total  radial  flow  in  velocity  space: 


r  Art  2  - 

G  =  Y  a 


where  a  =  eE/n  and  the  quantity  G  (the  gain)  represents  the  race  of  flow  of  electrons 
across  any  energy  value  in  the  direction  of  increasing  energy. 

7.3.1 

Taking  the  polar  axis  along  the  direction  of  grad  n  and  assuming  that  there  is 
rotational  syr.nacry  about  that  axis  so  that  the  associated  Legendre  polynomials  are 
not  required. 


The  Gradient  Term  ^ 


in  the  Boltzmann  Equation 
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in  w!i£ch  we  have  again  introduced  the  vector  notation  to  eliminate  the^  requirement 
on  the  polar  axis  P^. 


7.3.3  Magnetic  Terms 


Since 


div^(aij^  X  v)  =  0 


we  have 


div^(a}|^  X  vf)  =  (cuj^  x  v)  '  grcd^  f  =  v  V  (cuj^  x  grad^  f) 


The  expansion  inv'olves  the  associated  Legendre  functions  and  the  general  tern  will 
not  be  given.  The  first  two  terns  yield 


div  (co,  X  vf)  =  -V  •  (  '  X  f  ),'v 
V  b  n 


7 . 3.4  Expansion  In  Fourier  Series 


Ke  shall  assune  an  alterne.ting  electric  field. 


j.t  .  , 

e  ,  since  ve  can  alvays 


obtain  the  direct-current  formlae  by  setting  -.e  -  0  and  replacing  root-nean-squarc 
values  bv  the  dc  values.  li:e  distribution  function  vill  then,  in  general,  be  a 
Fourier  series  in  - 1 : 
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The  exponential  notation  Is  very  convenient  for  all  linear  equations.  Unfortunately, 

-*  I 

we  shall  have  to  take  products  such  as  a  f  and  the  meaning  of  ..such-a  pr-oduct  Is.  tlie_^ 
product  of  the  real  parts  and  not  the  real  part  of  the  product.  Performing  this 
operation  and  then  reintroducing  the  exponential  notation,  we  have 


2.  -  a  +  (2f‘  *  aJ  “  a  a»J 


I  I 

where  fj^^  Indicates  the  real  part  of  f^. 


7.3.5  The  Comnonent  Equations 


If  the  Boltzmann  integral  term  B  is  decomposed  as  well  as  the  electric  and 
magnetic  terms,  one  obtains  an  infinite  set  of  equations  of  which  the  zero-  and 
first-order  ones  are  given  below: 


V  , .  -t\  ,  _L_  d  /  2  -  vl  \ 

=  _  div  I  -r  ^  aT,  j 

3  0.2  d\,-  1  r 


0  “  3 


„0  .  -0  V  -rl  1  d  i  2  —  -rl  "Zl^  i 

=  JO  .  3  uiv  .  V  a  •  :  2 


V 


,  r,  dt, 

\  1  Ir 

1 

d  /  3  ,2  M 

/  ^  2  dv 

Sv^ 

dv  IrjJ 

4  =  j  a>  +  V  grad  +  5  +  a  |jd7v:^'o  ^  ^l) 

'1  d  3  /,-2  ^  ,2'\]  -  -*1 


In  the  case  of  electrons  (but  not  ions)  the  spherical  harmonic  expansion  con- 

2 

verges  well  and  the  terms  in  f  are  negligible.  The  Fourier  series  also  converges 

well  at  high  frequencies  (p^/v^  >  ^/2).  The  transition  to  low  frequencies  in  which 

(7-4) 

the  second-harmonic  terms  become  large  has  been  studied  by  Margenau  ;  this  case 
will  not  be  considered  here.  We,  therefore,  simplify  the  above  set  of  equations  to 


0  V  “^1  ]. 

®0  =  f 

6v 


“i  -  >  “  *  3  ^  ° 


■?1  1  r®  I 

Bq  =  V  grad  fo  +  2  "IT" 


=  j  "  fl  +  V  grad  f^  +  a  —  -  X  f^ 


The  divergence  of  the  alternating  flov  produces  an  alternating  electron  de-sity 
Oj^.  The  alternating  flow  is  dri%-en  by  a  n^  and,  hence,  n^  exists  in  virtes  of  a 
divergence  in  a  n^.  Since  a  is  divergeticeless,  n^  is  produced  by  the  interaction  of 
an  alternating  electric  field  and  a  concentration  gradient  in  the  same  direction. 
However,  if  such  situations  do  not  occur,  it  may  be  assumed  that  and  f”  are 
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.  .  r  1 


negllblble.  The  four  coaponent  equations  may  then  be  reduced  to  three  simpler 
equations: 


V  ,  Id 

®0=  3 

6v 


d  /  2  - 
2  dv  ( ^ 

V  N. 


Bq  =  V  •  grad  fg  -  X  f  Q 

-♦1  -♦  Q 

B-  =  J  05  £-  +  a  -  05.  X  f , 

1  •*  1  dv  b  I 


*  Q 
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Tha  Collisioa  Ip.coarals' 


7.4.1  Elastic  Collisions  Between  Electrons  and  Neutral  Particles 


7. 4.1. I  Geometry  of  a  Collision 


Consider  a  collision  betv/een  an  electron,  denoted  by  lower  case  letter,  and 


a  neutral  molecule,  denoted  by  capital  letters.  In  a  collision,  the  positions  of 


these  particles  do  not  change  appreciably  but  their  velocities  change;  we  shall 


3  3 

denote  by  primes  the  velocities  of  particles  about  to  be  scattered  into  d  v  d  V 


about  the  velocities  v  and  V.  We  do  not  consider  the  changing  velocities  during 


the  collision  but  only  the  velocities  before  and  after  the  particles  are  xvithin 
range  of  their  interaction  forces.  During  the  collision  the  velocity  of  the  center 


of  gravity,  v  ,  remains  fixed  and  if  the  collision  were  elastic,  the  relative  velocity 

s 


c  has  the  same  Mgnitude  before  and  after  but  has  changed  in  direction  by  an  angle 


X  ,  called  the  scattering  angle  in  the  center  of  gravity  syfetem.  The  center  of 

s 


gravity  divides  the  relative  velocity  inversely  as  the  masses,  so  that  the  velocities 


of  electron  and  molecule  remain  on  two  spheres  in  velocity  space  centered  on  the 


velocity  of  the  center  of  gravity  and  of  radii  Mc/(M+m)  and  mc/(M+m).  The  radius 


of  the  molecule's  sohere  is  very  small  compared  v.'ith  tli-at.of  the  electron,  so  that 


it  can  almost  be  considered  as  a  point  .  If  the  molecule '.ere  at  rest  before  the 


collision,  the  origin  of  velocity  coordinates  should  be  taken  at  V' 


(v  -  V')-  = 


r  n) 


-r  c  (I  -  cos  X  ) 

/  O’ 


or  if  V  =  0 


> 

v'  (I  -  cc 


r  n)‘ 


-  -  .  -  .,..{7-3)  ,  .  .(7-7) 

Parrs  or  ::>ecricn  3  arc  iron  Ai.113  ana  Oinaauri;  and  Curcvicn 


This  last  expression  gives  the  energy  loss  by  an  electron  in  striking  a  nolecale 
at  rest.  This  energy  goes  into  the  recoil  energy  of  the  molecule.  For  electrons 
it  is  sm-ill  and  we  can  write  it  in  the  form 


Av  Mn 

^  (M  +  m) 


r  (1  -  cos  X  ) 
2  g 


Another  useful  relation  from  the  triangle  V’,  v  ,  v  gives  the  scattering  angle 
in  the  laboratory  system  in  terms  of  the  scattering  angle  i**  the  center  of 
gravity  system.  It  is 


sln^  X, 


2  2 
sin  X 

_ _ g 

0  “  2  2 

M  +  m-  +  2M  m  cos  X 


The  scattering  function  oCX)  is  deduced  from  the  force  law  as  a  function  of 
^  3nd  c,  and  is  measured  experimentally  in  terms  of  X^  and  v',  and  it  is  always 
necessary  to  reduce  one  to  the  other.  The  dependence  of  cr  on  the  azimuthal  angle  ij' 
is  never  knotim  and  it  will  ba  assumed  that  a  does  not  depend  on  The  combination 
ccr  always  enters  the  equations  and  we  shall  use  the  sysbol  p(c)  =  cor  for  it  (p  has 
the  dimensions  of  a  volume  per  second  par  steradian). 


7. 4. 1.2  Reduction  of  the  Collision  Integral  Term 


The  rate  of  change  of  the  number  of  particles  in  a  velocity  volume  element  d  v 
is  given  by  taking  tlie  difference  between  the  number  of  particles  scattering  into 
that  volume  element  and  the  numher  of  scattered  out 


d^  V  -  I  F(V')  p  d-p.  dV  d\'  -  J  F(V 


(V)  f  (v)  p  d-p  d\ 


p  is  the  same  in  both  these  integrals  and  d  fi  is  the  solid  angle  of  the  scattering 


-.  3  3 

sphere.  F(V*)  d  V’  and  f(v')  d  v'  are  the  numbers  of  particles  in  ppsi.tion_tp,^ 


_  3  3 

scatter  into  the  volume  elements  d  V  and  d  v,  and  the  integrals  are  to  be  taken  over 


all  scattering  angles  and  all  molecular  velocities.  By  Liouville's  theorem 


3  3  3  3 

d  v'  d  V  =  d  V  d  V 


so -that  the  two  integrals  may  be  combined  into  one 


B  =  J  [F(V’)  f(v’)  -  F(V)  f(v)l  p  d^Q  d^V 


whfch  is  the  conventional  way  of  writing  Boltzmann's  collision  integral.  This  form 
is,  however,  rarely  convenient,  since  F(y*)  contains  V*  and  not  V  and  it  is  not 


easy  to  inegrate  it  with  respect  to  V.  It  is  more  convenient  to  keep  dV  and  convert 


3  3  3 

d  v*  into  d  v.  We  wish  to  determine  the  volume  element  d  v'  so  that,  for  given  V 


and  X  ,  the  final  velocity  after  the  collision  ends  up  in  d  v.  As  all  angles  are 
B  ' 


held  constant  and  V  is  held  fixed,  the  volume  element  d^v*  must  be  exactly  similar 
,3 


in  shape  to  d  v,  but  all  linear  dimensions  increased  in  the  ratio  {(v*  -  V')/(v  -  V')| 
so  that 


'  =  (jv'  -  V'i  [7  -  V'i)^  d^ 


This  loads  to  the  expression: 


/  F(V)  f(v*)  (i7’  -  V’lyrv  -  V'i)^  p  d-::  d\'  -  f  F(V)  f(v)  p  d'Q  d^V 
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7. 4. 1.3  Expnnsion  in  Spherical  Harmonics 


—r, 

We  now  substitute  for  f(v)  its  expansion  in  spherical  harmonics  ant!  for  f(v  ) 
the  corresponding  scries  in  terms  of  9'.  Using  the  addition  theorem  for  spherical 
harmonics,  we  have 


f(v*) 


P/O)  P,(X„)  +  2  J 


-^P>)  pX>  cos  S(9  -  (,) 


As  P  has  been  assumed  to  be  independent  of  the  azimuthal  angle  -Jr,  the  term  in 
cos  m(9-i|f)  goes  out  on  Integrating,  so  that  the  Boltzmann  integral  for  the  coef¬ 
ficient  of  P  (9)  becomes 


=  f  t'Cv')  f^V)  p^OCq)  p  d^n  d^v 


-  J  F(V)  f^(v)  p  d^n  d\ 


So  far  the  equations  have  been  perfectly  general,  but  we  must  now  make  assump¬ 
tions  which  in  practice  restrict  their  application  to  electrons,  or  at  most  to  light 
ions  in  a  heavy  gas.  We  shall  assume  that  the  mass  ratio  m/H  is  zero,  so  that  v'  -  v 
and  that  F(V)/ng  is  a  5-function.  Then, 

y  =  0  =  0 

y  0  b"  =  -f‘(v)n  r  p(v,X„)[l  -  P.(X-)]  d^a  =  -V  ,f^ 

g  J  U  ^  U  Cr. 


96 


7.4. 1. 3  E::p;ir!sion  in  Sph.;rical  Har.-noaics 

We  now  sjbsCitute  for  f(v)  its  expansion  in  spherical  harmonics  and  for  f(v*) 
the  corresponding  series  in  terms  of  3'.  iJsing  the  addition  theorem  for  spherical 
harmonics,  we  have 


f(v') 


■  E ''  h 


(9)  P.(x„>  *  2  ^  P^(0)  pjex^)  cos  Sw  -  t) 


As  p  has  been  assumed  to  be  independent  of  the  azimuthal  angle  Jr,  the  term  ^n 
cos  goes  oat  on  integrating,  so  that  the  Boltzimnn  integral  for  the  coef¬ 

ficient  of  P  (9)  becomes 


-  J  F(V)  f^(v)  p  d^fl  d\ 


^  So  far  the  equations  have  been  perfectly  general,  but  tje  iqust  how  make  assumo- 

trons  which  in  practice  restrict  their  application  to  electrons,  or  at  most  to  light 

ions  in  a  heavy  gas.  We  shall  assume  that  the  mass  ratio  m/M  is  zero,  so  that  v*  «  v, 

and  that  F(V)/n  is  a  5-functioa.  Then, 

g  ’  - 

d  =  0  3*^  =  0 

i  4  Q  ■  =  -f‘(v)ng  /  p(v,Xq)[i  -  p  .(Xq)]  d^a  =  -v^;f" 


Tp.e  "collision  freqaencias"  ,  introduced  above  are  defined  in  the  laboratory 
reference  systeni.  Introducing  an  arbitrary  definition  the  first  three  collision 

frequencies  are 


n  p  da 
g 


n  p(l  -  cos  X^)  d  a 


3  2  2 

“g  p  i  sin  Xq  d  a 


The  total  collision  frequency  is  v^q.  The  integral  is  generally  improper  as 
p  has  a  singularity  in  the  forward  direction.  The  apertures  of  experimental  ap¬ 
paratus  prevent  observations  from  covering  this  singularity,  so  that  the  quantities 
tabulated  as  experimental  "probabilities  of  collision"  are  expressed  in  terms  of  an 
incomplete  integral. 


The  collision  frequency  for  "momentum  transfer",  is  of  great  Importance 

in  transport  theory  and  the  simpler  symbol  is  used  for  it.  It  is  a  weighted 
collision  frequency  in  Which  the  backward  scattering  counts  double,  right  angle 
scattering  has  the  weight  one,  and  forward  scattering  is  not  counted.  The  frequency 
might  well  be  called  the  "scattering”  frequency;  here  right-angle  scattering  has 
the  weight  3/2  and  neither  forward  nor  backward  scattering  is  counted. 


7. a. 1.4  Recoil  of  the  Heavy  Molecule 

The  conclusion  that  =  0  means  that  collisions  do  not  change  the  energy 
distribution,  and  it  follows  irem  the  assumption  that  the  molecules  were  clastic 
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O 


I 


li  =  (v  V  r) 

d(v2) 


Hejcc  we  take  account  of  the  Mixwell  distribution  of  the  molecules  F(V'),  v;hich 


was  previously  replaced  by  a  delta  function,  and  this  will  be  done  in  a  somewhat 


qualitative  way.  In  integrating  over  d  V*  the  point  V  moves  over  a  MaxiJell  distri¬ 


bution  centered  about  V'  =  0.  For  fixed  v  and  scattering  angle  X  ,  this  means  that 

£ 


v'  moves  about  v^,  its  value  for  V*  =  0,  by  amounts  of  the  order  V*. 


v'  =  7^ 


f(7')  =  f(v')  +/(v'^)  — ^ 

a(v'2) 


when  tfiis  correction  is  introduced  into  the  expression  for  the  collision  integral. 


one  obtains 


„0  2m  I  d  i  3  /'-O  .  2.  df''  \ 


df®  Vi 


But  >..'e  k Lit)-.;  that  if  the  electrons  have  a  Maxwell  distribution  at  the  same  temperature 


as  that  of  the  g.as,  collisions  v.ill  not  change  t!i  is  distribution  and,  therefore,  this 


last  eauatioa  mast  give  zero  if  f  is  a  M-ixi:ell  distribution.  This  determines  the 


factor  to  be  equal  to  2k  T'_/iu,  and  we  have 


.0  2lv! 


2k  T  .  .0  - 

g  d  l  i 


for  Eho  rate  of  chrmre  of  Ciie  energy  distribution  of  electrons  of  m-ass  m  owing  to 
clastic  collisions  with  molecules  of  finite  miss  M  having  a  temperature  T  . 


7-4.2  Elastic  Collisions  of  Electrons  with  Ions 


To  describe  elastic  collisions  between  electrons  and  ions  one  can  employ  the 
foregoing  general  expressions  for  the  integral  of  elastic  collisions  betv;aen  an 
electron  and  neutral  particles,  vjithout .modification,  since  die  only  assumption  made 
in  their  derivation  vns  m  «  M.  It  is  necessary  only  to  calculate  the* number  of 
collisions  between  the  electrons  and  the  ions,  y^(v).  For  this  purpose  one  sub¬ 
stitutes  in  the  collision  iuLegral  the  Rutherford  formula  for  the  differential  effective 
cross  section  for  the  scattering  of  an  electron  by  an  ion.  We  then  have 


v.(v)  =  2nN  v(c^/2.nv^)^  f  - — sin  9  dS 
1  1  I  .  H  O 

q\  2 

min 


=  2nN.  ■  -u  £n  f  I  +  cot 
V  2  3  V 

m  V  V 


where  N.  is  the  concentration  of  the  ions,  which  for  simplicity  are  assumed  to  be 
sinulv-char-ued . 


If  '..’c  consider  the  scattering  of  an  electron  by  a  free  ion,  then  integration 
s'ncald  be  carried  out  free-.  0  to  ::  (i.e.,  J  .  --  0)  and  the  collision  frequency 

diverges  Logarici'.nical  Ly  at  sc-a  1 1  0 .  In  a  plascia,  luivever,  the  ions  are  not 
entirely  free:  as  a  r-esait  of  interaction  beL'..een  the  ions  and  the  electrons,  the 
field  of  e.ic’.i  ion,  und.er  equ  i  I  ibr  uir:i  conditions,  has  a  Coulomb  character  only  to 


on  the  I'rder  of  tb..-  beb ve  radius  D: 


O 


r  kTkT  ni 

D  .  ^ -  1 

■4ne  N(kT  +  kT  )  J 


where  T  =  ion  temperature  and 

t^  =  electron  temperature. 

At  distances  greater  than  D,  the  Coulomb  field  of  the  ion  drops  off 
as  a  result  of  screening.  Consequently,  D  is  the  laaximum  distance  at  which  a 
substantial  interaction  between  the  electron  and  the  ions  still  takes  place,  i.e., 
the  maximum  Impact  parameter.  It  can  be  used  to  express  the  minimum  scattering 
angle: 


d  2  tan‘^(e^/mv^D)  S:  2e^/mv^D 

min 


Therefore: 


,  X  -  e^  . 

V  (v)  =  2nN  —  in  U  ^ 

m  V  b  e  / 


2  2  ^  3  3  6  X 

It  is  important  that  D  (kT^)  e  ~  k  e  N  al\«y3  be  a  large  quantity  in 

the  cases  of  interest  to  us.  This  means  that  the  second  term  in  the  logarithm  is 
alv.’-iys  the  principal  term.  Consequently,  the  principal  contribution  to  the  number 
of  collisions  betv/een  the  electron  and  'the  ions  is  made  by  the  weak  scattering  — 
scattering  by  small  angles.  In  one  such  collision,  the  change  in  either  the  energy 
or  in  the  electron  momentum  is  insignificant.  In  fact,  the  fraction  of  the  energy 
lost  by  the  electron  v.hor.  scattered  by  an  angle  3  is  5^  -  2ri  (1  -  cos  3)/M.  Consid¬ 
ering  that  the  principal  role  is  played  by  collisions  chnt  lead  to  the  scattering 
by  a  small  angle,  on  the  order  of  0  .  , 


we  fixid 


,  /  2,.1/J  3 

O,  6,  =  t#  .  w  =  M  I  f.T  j  « 

k  k  .  Bin  M  M  \  kX  i 

Bin  e  ' 


Analogously,  the  change  in  the  noasutum  Is 


6^.  Ip-Pll  / 


It  mast  be  emphasized  that  although  the  change  in  the  momentum  In  one  Impact  Is  small. 


the  change  in  the  energy  is  considerably  smaller:  ~  m/M. 

K  P 


Scattering  by  large  angles  adds  to  the  number  of  collisions  only  a  term  of 
order  unity,  which  is  small  compared  with  the  main  logarithmic  term.  A  similar 
correction  in  equilibrium  plasma  results  from  an  exact  solution  of  the  problem  of 
scattering  in  a  Debye  field. 


Inelastic  collisions  between  electrons  and  ions,  which  lead  to  their  excita¬ 
tion  and  multiple  ionization,  do  not  differ  at  all  from  inelastic  collisions  with 
neutral  particles,  considered  above.  However,  oxving  to  the  large  values  of  the 
maximum  elastic  impact  parameter  (2),  the  role  of  inelastic  collisions  is  greatly 
reduced.  Collisions  of  electrons  accompanied  by  bremsstrahluag,  which  is  of  im¬ 
portance  at  high  electron  energies,  wilL  not  be  considered  here. 


7.4.3  Electron-El-jclron  Collisions 

T'ne  principal  role  in  a  collision  between  aa  electroa  aad  ions,  as  seea  above, 
is  played  by  long-raage  collisions,  which  lead  to  weak  scattering.  Both  the  energy  an 
the  momentum  of  the  electrons  are  changed  only  slightly  by  o.ne  such  collision.  This 
is  the  conseqs’.ence  of  the  -singularity  of  do  Coulomb  interaction  and  therefore  pertains 


not  only  to  collisions  between  electrons  and  ions,  but  also  to  collisions  between 
electrons  and  electrons.  The  difference  lies  only  in  the  fact  that  the  fraction 
of  the  energy  and  the  fraction  of  the  momentum  lost  by  the  electron  when  colliding 
with  another  electron  is  of  the  same  order  ~  1*  whereas  in  the  collision  with 

the  ion  5.  75  ~  m/M. 

*  P  ■ 

Thus,  in.  considering  the  integral  of  inter -electron  collisions,  one  can  use  the 
differential  expressions  derived  earlier  for  B.  In  addition,  we  can  integrate  this 
expression  over  the  scattering  angles  dfl  (using  the  fact  that  o'(0,u)  has  a  sharp 
maximum  at  5  ~  0).  We  then  find  that  the  integral  of  the  collisions  between  the 
electrons  are  satisfied  by  the  expression 


B  =  -  div  j 
ee  V 


where 


-Jy  ^  ^  f  dvj^v(u)  {u^(f(v)  grad  ^  f(v,)  -  f(vj^)  grad  ^  f(v)) 


-  u[f(v)  (u  grad  ^  f(Vj^))  -  f(Vj^)(u  grad  ^  f(v)]) 


and  is  the  velocity  of  the  particle  with,  which  the  electron  collides  and 
u  -|v  -  '.‘(u)  is  the  number  of  collisions,  where  v  must  be  replaced  by  u  and 

IJ  by  account  is  also  taken  of  the  fact  that  the  scattering  particles  are 

electrons,  i.e.,  that  F(vj^)  ^  f(Vj^). 

Let  us  consider  c.oj  B  -  the  integral  of  intcr-e lectron  col  I  i« ions  for  the 

ee 

function  f  .  by  putting  f  =  into  t!ie  expression  for  the  collision  integral 

and  ir.teorutiug  over  the  angles,  one  obtains: 


“cc '  -  TS;  "o  *  w 

V 


where 


“  i  f  ^1  ^  ‘  cos^  ep  v(u) 


=  f  vj  fo(7^)  dv^ 

0 

*2  =  J  J  IfJ  f|)  (1  -  ‘o«^  v(u)  dVj 

-  /  '’}  ^0<’'l)  “’'l  *  ’'"  /  ’l  'O*'!’  ‘"'l 


Hare  9j^  is  the  angle  between  v  and  Vj^,  u  =  |v  -  v^l;  when  integrating  over  the  angles 
we  neglected  the  variation  of  the  logarithmic  term  in  v(u)  "(compared  with  the  varia- 
tion  of  the  principal  term  ~  1/u  ).  For  fast  electrons,  whose  velocity  v  is  much 
greater  chan  the  average  velocity  of  the  plasma  electrons,  the  coefficients  Aj^  and 
assume  a  simple  form:  Aj^  =  v(v),  and  A^  =  2K  v(v)/3'n,  where  K  is  the  average  energy 
of  the  scattered  electrons  (in  Che  case  of  a  Maiwellian  electron  velocity  distribution, 
2K  /3a  -  kT  /m). 


.4 . 4  Inelastic  Electron-Neutral  Collisions 


Inelastic 
Che  excitation 
In  addition,  s 
excited  state 


collisions  between  electrons  and  neutral  particles  are  accompanied  by 
of  rotational,  vibrational,  or  electronic  levels,  and  also  by  ionization, 
o  called  second-order  impacts  are  possible,  in  which  the  energy  of  the 
of  the  molecule  is  transferred  to  the  Inco.iing  electron.  .4n  exact  cnl- 
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culaclon  of  all  these  inelastic  processes  is  quite  complicated;  in  addition,  their 
cross  sections  are  kno-Am  only  in  a  few  cases.  Therefore  there  is  no  complete  theory 
of  inelastic  collisions  In  which  the  problem  is  solved  as  accurately  as  in  the  case 
of  elastic  collisions.  In  spite  of  this  fact,  it  is  possible  to  analyze  relatively 
simply  two  important  limiting  cases;  specifically,  we  consider  cases  when  the  electron 
energy  is  considerably  greater  than  the  energy  of  the  e:ccited  level,  or  the  ionization 
energy  (K  » "k  o),  and  when,  to  the  contrary,  the  electron  energy  is  only  slightly 
higher  than  the  excitation  energy  (K  -  iko  «  Kand  positive). 

In  the  former  case  the  expression  for  the  integral  of  inelastic  collisions  is 
found  in  the  same  manner  as  for  elastic  collisions.  It  is  merely  necessary  to  con¬ 
sider  that  the  energy  lost  by  the  electron  in  inelastic  impact  is  consumed  essentially 
in  excitation  of  the  molecule,  and  this  is  connected  with  a  transfer  of  energy  "feu 
(thus,  v*  -  V  =  hco/mv);  in  this  case  the  neutral  particle  simply  goes  from  the  ground 
state  into  the  excited  state.  We  then  have 

1  O  2  _ CD  _ 0  ,  c 

inelastic  ~  _  2  ^v  ^  ^’co  m  dv  0 

2v 

^inelastic  ^1 

Here  v  is  the  number  of  inelastic  collisions  accompanied  by  the  excitation  of  a 
quantum  ki)  (as  we  shall  call  the  transfer  to  the  molecule  of  an  energy  ka,  consumed 
in  excitation  of  soma  level)  ^ 

(v)  =  v(N^  I-  (v,  3)  (1  -  cos  3)  du 

io  •  m  m  /  CD 
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where  is  the  differential  effective  scattering  cross  section  in  inelastic 

0 

collision.  N  and  N  are  the  nu.uber  of  molecules  in  the  ground  and  excited  states 
’am 

respectively.  Furthermore,  r^Cv)  is  the  fraction  of  the  energy  lost  per  unit  time 
by  the  electron  to  excitaticn  of  a  quantum  to 

r  (v)  =  -  N®*  )  V  f  cr  (v,d)  dQ 

fio  2  m  .a  /  £0  ’ 

fflV  — 


and  T  is  the  effective  ternerature 
to 


.  N°  +  N®’' 
X  _  SS  m  m 
to  ~  2k  „0 


N  +  N 
ffl  m 


ex 


It  is  Important  to  emphasize  chat  in  the  case  when  the  quantum  to  is  small  not  only 
compared  with  the  electron  energy,  but  also  compared  with  the  energy  of  neutral 
particles  (Ko  «  kX),  and  if  the  neutral  particles  have  a  Boltzmann  distribution 


N®^  /N?  =  expC-iko/kT)  (i.a.,  if  Che  collisions  with  the  electrons  do  not  substantially 
mm 


change  the  number  of  excited  molecules),  then  the  effective  temperature  is  equal 


to  the  molecule  temperature  Xj^. 


In  the  second  limiting  case,  when  the  electron  energy  is  slightly  greater  than 
the  excitation  energy,  the  coILiJing  electron  merely  goes  from  the  region  of  large 
energies  into  the  region  of  snail  energies  (X  0).  Therefore  at  large  energies 


where  v  is  the  total  frcqjency  of  the  excitations  of  the  level  (it  is  assumed 

that  «  N*^).  Tlie  fact  that  the  electron  cannot  simply  vanish  but  goes  into  the 

mm 

region  of  small  energies  (K  ~  0)  Is  taken  into  account  here  by  adding  to  the  equation 

2 

for  £q  a  6-functlon  source  of  electiois^  -Q6(0)/4nv  , 

00 

where  Q  =  dN/dt  =  4it  F  f^  dv. 


We  note  in  conclusion  that  in  the  general  case  the  characteristic^ dependence 

of  the  total  cross  section  of  the  inelastic  collision  on  the  electron  energy  is 

such  that  the  total  cross  section  aCv)  vanishes  when  K  < 'feoj,  then  it  increases, 

reaching  a  maximum  X  ~  (3  to  5)  iioi,  and  then  starts  to  diminish  slowly.  Since  in 

the  average  impact  K  » "W  and  the  electron  loses  an  energy 'tra,  the  electron  in  the 
most  probable  inelastic  collisions  loses  a  small  fraction  of  its  energy. 

We  note  also  that  in  those  cases  when  not  one  but  several  levels  iia)  can  be 

excited,  we  have  B  =  ZB  , ,  It  must  be  taken  into  account  also  that  some  "Of  the  in- 

,  oi 

X  • 

elastic  collisions  --  lonixatlon  and  the  effective  recombination  (recombination, 
capture  of  electron  by  a  molecule,  etc)  —  are  accompanied  by  a  change  in  the  number 
of  electrons  in  the  plasma.  It  is  therefore  necessary  to  add  in  the  collision 
integral  for  the  function  f^  the  terms 


VA./ 

-Vj(0  fu*  f  £„(v')  <lv’ 

V  m 


the  first  of  v;hich  describes  the  effective  recombination,  and  the  second  describes 

ionic, iCion.  Here  (v)  is  the  total  recombination  frequenev 

r 


and  V,  (v 
ion 

Is  the  Ionization  frequency,  i.e.,  the  number  of  ionizations  produced  per  second  by 
electrons  of  velocity  v',  which  lead  to  the  appearance  of  a  new  electron  of  velocity  v, 
where  is  the  ionization  energy.  These  terms  usually  do  not  exert  a  noticeable 
influence  on  the  form  of  the  distribution  function;  they  do  determine,  however,  the 
concentration  of  the  electrons  in  the  plasma. 


’,v) 


=  N  v*  fa. 
m  J  ion 


(v‘,v,i-)  da 


lOS 


-Solutions  of  the  Boltzmann  Eauat: 


ion 


If  a  spherical  ha^aonic  expansion  of  the  electron  distribution 


function  is 


/  v.fi. 

^  ~v  )*  Boltzaiaa  equations  for  the  zero  and  first 

order  distribution  functions  (taking  into  account  electron-electron,  elastic  electron 

are:  • 


ion  and  elastic  and  inelastic  electron-neutral  collisions) 


Ms 

at 


en 


-  .(v)fj  -  (ep 


Here, 

‘eo'V  ■  *  7  f  Aj(y  ^  J 


,  4iti!(v)  2.  .  . 

\  =  ~N  J  ^1 


.  Anv(v)  f  2  .  ^  "1 


EL 


-  .2s 

M 

-  e'f 


^ on-neu^*. al  olastLc  cc._iSLon  frequency 
1^— '-^on— ion  eluscic  coliisxon  frequency 


Ei.VS 


..i:.'EL.\S 


V  is  Che  total  nuabcr  of  cluctrj:n  aoliisioas,  vhere 
-  -55;™  )  j\.  Cy,-)(l-cosg)dS 


(fo)  is  the  collision  integral  for  the  function  f^,  which  describes  in¬ 
elastic  collisions  between  electrons  and  molecules. 

It  should  be  brought  out  that  the  electron  energy  relaxation  time  is  always 
much  greater  than  the  momentum  relaxation  time.  This  implies  that  the  relaxa¬ 
tion  time  is  much  greater  for  the  function  f©  than  for  the  function  fj.  As  a 
consequence,  the  function  f^  alv/ays  changes  more  slowly  with  time  than  the  func¬ 
tion  fj.  Hence,  when  the  equation  for  fj^  is  integrated  with  respect  to  time, 
the  function  fj,  can  be  considered  time  independent  to  a  first  approximation. 

In  the  equationfbr  the  function  fg,  the  last  term  on  the  right  side  of  the 
equation  (®gg(^o)  )»  which  represents  electron-electron  collisions,  is  of  the 
order  u  f  ,  where  is  the  4iter -electron  collision  frequency.  The  remaining 

c  O '  ® 

terms,  which  describe  the  collisions  between  electrons  and  heavy  particles,  is  of 
the  order  62  f^,  where  6  =  ^™/m,  u  =  Hence,  the  relation 

between  v  and  5.?  determine  the  form  of  the  function  f^.  A  highly  ionized  plasma 
is  characterized  by  the  condition  y^>>6u;  and  a  v/eakly  ionized  plasma  by  the 
condition  «Sv(in  a  fully  ionized  plasma', '  .  ’ 


7.5.1  Hi.^’nlv  Ionized  Plasma 


In  a  ; 
tion  funzt 


cral  Eoltr- 


colliaic-'.s 


00  > 


.ghlv  ionized  plasma  (v^  »  5v) ,  the  form  of  the  zero  order  distribu- 
n  f j  is  set  by  the  inter-alectron  collisions.  The  solution  of  the  gen- 
equations  may  be  obtained  by  the  method  of  successive  approximations 
-^...  where  the  zero  order  approximation  includes  only  inter-electron 
If  the  Boltzmann  equations  are  solved  for  the  zero  order  approximation 
aini  a  Max\.‘ellian  distribution  function  in  the  case  of  a  hem.ogoneour. 
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plasm: 


f  =  /— S_j 


"•oo  ‘•v2nkT 


exp 


f J 


The  Maxivellian  distribution  becomes  established  within  a  time  Vu-,.  When 


V  »  cu  this  process  is  much  faster  than  the  process  of  transfer  of  energy  to 


heavy  particles,  this  implies  that  f  is  close  to  being  Maxwellian. 


If  the  expression  f^  =  £qq  +  +...  is  substituted  into  the  Bolbzmanh 


equation  for  f^,  and  then  both  sides  are  multiplied  by  4ttv  and  integrated  with 
respect  to  v,  one  obtains: 


3?  r  *  f  4-^  .  0 


If  the  Max^Jellian  distribution  is  taken  for  f^Q,  then 


dN 


-rr^  +  [w  -  u.  ]  N  =  0 
dt  rec  ion  e 


where 


>  '  ... 

i..  =\/Tf  r./S*  v^v')2,  (v,v-)  exo  dv- 

lon  „  V  kX^  y  Jo  V  m  im  ’  *  I  2kT^  j 


and 


m  'N 


3/2 


rec 


kx  y 

e 


)  dv 


7,5.2  W-jaklv  Ionised  Plasm 


In  c'-.e  case  of  a  '.eiikly  ioniced  nla.s.ru,  the  collisions  bef-cen  elect 


in.siyn  It  leant  in  tne 


;it2f-inn  equation  Cor  f  ^  (since  v  ,  and  those 
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sions  may  be  disregarded  in  a  first  approximation.  laey  are  even  less  signifi¬ 
cant  in  the  equation  for  fj^,  since  «  by  «  y.  Hence,  the  function  fj  in  a 
homogeneous  weakly  ionized  plasma  (to  terms  of  order  5)  is  given  by  the  expression 


f  _  _ o 

^1  -  "  “  ^ 


where  u  satisfies  the  equation: 


ou  ,  /  .  eE  .  e  r-^  T 

+  i;(u)  u  =  —  +  —  [uxHj 
dt  m  me  o-* 


u  is  the  velocity  of  the  directed  motion  of  the  electron.  If  this  expression  for 

fj  is  substituted  into  the  general  Boltzmann  equation,  the  following  equation  for 

f  is  obtained; 
o 


^o  1  3  C  ^  [  ('k  j.  \ 

‘v  +  "i'  •; 


kT  2eEu 
m  '3m 


N 

5  v( 

y  ov 


IHELAS  ^  p 

ne  '  O' 


Depending  on  the  relation  between  the  tine  1/x,  dirring  which  the  electric 
field  changes  significantly,  and  the  relaxation  tine  fer  the  function  ~  l/£a 

we  distinguish  here  cases  of  slov/ly  varying  field  (jc  «  :•/,  and  rapid  ones 
(ju  »  5i;)  (the  sane  as  in  the  amilysls  of  the  electron  temperature  in  elementary 
theory  or  in  a  strongly- ionized  plasr-j).  In  the  former  case,  which  is  quasi- 
statiunary,  the  dependence  of  f  on  the  tim.e  can  be  disregardevi ;  in  particular. 


this  takes  place  r-'.curally  in  the  case  of  a  eoastjiiL 


field.  On  the  other 


hand,  for  a  rapidly  alternating  eK:Ccric  rield,  >■  »  t  the  function  doee;  not 


have  a  chance  to  change  as  rapidly  as  the  field;  it  therefore  settles  at  a  certain 
average  level,  independent  of  the  cisie,  and  the  variable  deviations  from  this 
level  are  small,  of  amplitude  on  the  order  of  6y/x  (the  same  as  the  observations 

of  the  electron  temperature  in  elementary  theory).  Consequently,  in  both  cases 

# 

we  can  neglect  in  first  approximation  the  term  df^/St,  and  thereby  get  rid  in 
fact  of  the  time  variable.  This  allows  us  to  find  an  analytical  solution  for  the 
above  equation  for  many  Important  cases:  for  elastic  collisions, _ in  inert  gases 
and  in  a  molecular  plasma.  We  now  proceed  to  analyze  these  solutions. 

7.5. 2.1  Elastic  Collisions 

If  all  the  collisions  are  elastic,  =  0.  Therefore  in  a  constant 

ne 

electric  field  E  we  have,  u  =  eE/mu  at  H  =0  and  the  above  equation  becomes; 

o 


+  6  wf^ 


1] 


1  a 
7^ 


(v^J  )  =  0 


1  2 
(Hera  v  =  I7(v)  =  .  Multiplying  this  equation  by  v  and  integrating  from 

2 

0  to  V,  we  see  that  =  0,  since  in  the  absence  of  an  electron  source  [v  0. 

Integrating  now  the  equation  =  0  over  the  velocities,  we  obtain 


f 

o 


C 


V 


mv  dv 


2  2 
2e  E 

3m5  'j" 


We  therefore  obtain  a  Maxwellian  d 


field  the  distribution  function  f 


istributlcn 
may  differ 


in  a  weak  field,  but  in  a  strong 
substantially  from  Maxwellian,  since 


o 


j  depends  on  v.  For  e:-:a:nple,  in  i  strong  electric  field  upon  collision  with 
molecules  -  hard  spheres  -  the  function  f^  is  determined  by  the  well  known 
Druyvesteyn  formula 


2 

where  £  =  v/i;(v)  =  l/na  is  the  mean  free  path  of  the  electron,  C  a  constant 
determined  from  the  normalization  condition,  and  the  term  kT  is  neglected,  this 
is  permissible  for  a  strong  field. 

The  Druyvesteyn  distribution  at  large  electron  velocities  differs  greatly 
from  Maxwellian;  it  drops  off  much  more  rapidly  than  a  Maxwellian  one.  The  cal¬ 
culation  of  the  function  f^  with  allowance  for  the  exact  dependence  of  the  colli¬ 
sion  frequency  and  the  velocity  for  'different  inert  gases  was  made  in  references 
7-10  and  7-11.  The  effect  of  a  constant  magnetic  field  is  taken  into  account  in 
reference  7-12  [the  magnetic  field  changes  the  velocity  v(u),  and  accordingly 
f^  also  changes]. 

Wa  considered  above  only  the  case  of  a  constant  quasi-stationary  electric 
field  ( t  »  cv)  icr  in  tr.is  case  we  can  neglect  in  first  approxi¬ 
mation  the  derivative  of  /jt.  The  function  f  now  assumes  the  forra 

o  o 


( 

fo  =  c  exp  -  , 

C  O 


mv  d\ 


*  V 


)  -  (v) 


1 


the  function  v(v)  without  the  r.ugnetic  field  is  equal  to  [.s' 


-1 


a  lid  i  ri 


lijid 


2  2  2 
»  cos  0  .  sin  0  ,  sin  3 

?(v)  *  -5 — — ~  + - 5 - 5-  + - 5 - 5- 

2[(a-ti^)  +  u  ]  2[(w-uUjj)^  .+  v  ] 

where  0  Is  the  angle  between  £  and  H,  oUjj  Is  the  gjrromagnetic  frequency,  is  the 
amplitude  and  ui  is  the  frequency  of  the  alternating  electric  field.  Such  e;lastlc 
collisions  are  produced  in  monatomic  gases  at  low  electron  energies  (up  to  1  ev). 
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1. 5,2.1  Inelastic  Collisions  in  a  Molecular  Plasria 

A  molecular  plasma  may  fae  defined  as  one  formed  by  diatomic  or  polyatomic  gases. 
In  such  a  plasma,  not  onfy  ekctrcnic  levels,  but  also  rotational  and  vibrational  levels 
can  be  excited,  the  energy  of  such  levels  being  exceedingly  low  (jicu  ~  10  ^  to  10  ^  ev 
for  rotational  levels,  and  iku  %  0.1  to  0.5  ev  for  vibrational  levels).  Hence,  in¬ 
elastic  collisions  in  such  a  plasma  become  important  at  electron  energies  of  the 
-2  ■ 

order  10  ev  (room  temperature). 

For  a  plasma  composed  of  a  diatomic  gas  such  as  hydrogen,  oxygen,  or  nitrogen 
at  electron  energies  less  than  1  ev, the  principal  role  is  played  by  losses  due  to 
the  excitation  of  viTrntional  levels.  The  energy  of  these  vibrational  levels  is  small 
compared  with  the  average  electron  energy.  For  such  cases,  the  integral  for  inelastic 


collisions  for  the  function  f^  can  be  represented  in  the  form: 


.inelas 


/e  \  1  ^  i  r  X  fkT  ^^0 


+  V  f , 


where  Ry(v)  =  T.  r^  describes  the  energy  losses  of  electrons  in  inelastic  collisions. 


When  this  expression  for  B 


Inelas 


is  substituted  into  the  Boltzmann  equation  for  fg 


for  a  weakly  ionized  plasma  (v^  «  5v),  it  may  be  shown  that  the  resulting  equation 
coincides  with  the  equation  considered  above  for  the  case  of  elastic  electron-neutral 
collisions.  It  is  merely  necessary  to  make  the  transformation: 


V,)  + 


elas  m  inelas 


For  the  case  of  a  strong  constant  electric  field,  instead  of  a  Druynesteyn  distribu¬ 
tion  which  one  obtains  for  the  case  of  elastic  electron-nCutral  collisions,  the  zero 


order  distribution  function  for  a  molecular  plas-ia  is  given  by: 


fg  *  c  exp 


In  order  to  obtain  an  explicit  expression  for  £q,  it  is  necessary  to  calculate  the 
function 


el 


6(v)  = 


5gi(v  +  V^)  +  Rjj(v) 
v(v) 


This  may  be  found  If  the  function  Ru(v)  =  Z  r  is  known.  However,  to  find  r  it 

^  i  “i  “i 

is  necessary  to  know  all  the  cross  sections  of  the  various  inelastic  processes. 

Seaton^^  has  calculated  the  cross  sections  for  the  excitation  of  electronic 
levels  by  electron  impact  for  atomic  oxygen  and  atomic  nitrogen,  Meyerott^^”^^^  has 
given  the  cross  sections  for  excitation  by  electron  impact  of  some  of  the  levels 
of  molecular  nitrogen.  The  exchange  of  energy  between’ electrons  and  the  rotational 
motion  of  N2  has  been  discussed  by  Gerjuoy  and  Stein^^  The  interaction  of  the 

electron  with  the  quadrupole  momant  of  the  molecules  was  considered.  Closed  expres¬ 
sions  for  cross  sections  for  an  electron  producing  a  charge  of  rotational  quantum 
number  AJ  =  ±  2  as  functions  of  J  and  of  the  ratio  of  initial  and  final  electron 
momanta  are  given.  If  it  is  assumed  that  the  average  electron  energy  is  much 
greater  than  kT^  (where  is  the  gas  temperature),  the  fractional  energy  loss  per 
collision  may  be  regarded  as  snail  compared  to  unity  for  most  electrons.  Under  this 
assumption  and  assuming  (the  probability  of  the  mvalecale  being  in  rotational 
level  J)  is  small  when  J  is  sm-ill,  the  cross  section  becomas: 


•>  9 

4a  Q-  a- 


where  Q  =  molecular  quadrapole  moment  In  units  of  rt 

Q  =  0.96  for  nitrogen 
ajj  =  first  Bohr  radius 

Haas^^  and  Schulz^ have  shovm  that  there  is  a  strong  probability  for 

electrons  of  1  to  4  ev  energy  to  lose  energy  equal  to  one  or  more  vibrational  quanta 
in  These  authors  interpret  the  process  as  the  formation  of  an  unstable  ion 

which  can  decay  into  various  vibrational  levels  of  N2.  Haas'  experiment  was  con¬ 
ducted  with  a  swarm  tube  and  demonstrated  that  the  maximum  of  the  cross  section  for 

•16  2 

this  energy  loss  process  occurs  at  2.3  ev  and  is  about  3  x -10  cm  in  magnitude. 


There  is  some  evidence  that  electrons  also  excite  vibration  in  0^  at  energies  of  a 
few  tenths  of  an  electron  volt.  Because  the  investigation  of  this  process  by 
mobility  and  diffusion  processes  is  complicated  by  attachment,  it  is  difficult  to 
assign  cross  sections  to  specific  processes. 

Since  the  cross  sections  for  the  various  inelastic  processes  that  may  occur  in  a 

(  7-7) 

molecular  plasma  are  not  all  well  kno;vn,  Ginzburg  has  adopted  another  method  to 

determine  the  total  loss  function 


R(v)  =  +  Rjj(v)  =  v(v) 


The  fraction  of  the  energy  lost  by  the  electron  in  a  highly  ionized  molecular  plasma 
is  related  to  cite  function  RCv): 
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The  above  relation  may  be  considered  as  an  integral  equation  for  R(v),  since  the 
product  ®eff  known  from  experiment.  Once  R(n)  is  determined,  6(v)  =  R(v)/v(v) 

is  known.  The  results  of  such  a  calculation  for  hydrogen,  oxygen,  nitrogen  and  air 
are  given  in  Figure  7 -1. 

If  the  calculations  for  6(v)  are  substituted  into  the  expression 


=  c  exp 


1-/- 


the  distribution  function  for  electrons  in  a  molecular  plasma  is  determined. 

The  results  of  this  calculation  for  electrons  in  hydrogen  in  a  high  frequency 

2  2  2 

electric  field  are  sho-.m  in  Figure  7-2,  The  ordinates  represent  -v  /v  where  v  =  2Ic'l7n; 
is  the  mean  square  electron  velocity.  The  dotted  line  represents  a  Maw^ellian  dis- 
tributj.on  (the  distribution  would  be  Maxijrellian  if  0  were  Independent  of  v,  as  it  is 
for  elastic  collisions).  It  mny  be  noted  from  Figure  7-2  that  for  electrons  in 
the  deviations  of  the  distribution  function  froiaMaxwellian  are  not  large, 

(  7**  18) 

Recently,  Caldirola''  has  determined  the  velocity  distribution  function  for 

electrons  in  a  spatially  homogeneous  slightly  ionized  plas-n  under  the  action  of  a 
constant  nacnetic  field  (i,e.,thacof  the  Earth)  and  an  alternating  electric  field 


O 


€ 


(the  tine  varying  niagnetic  contribution  to  the  force  exerted  on  an  electron  is 


assumed  negligible  compared  with  the  electric  one).  Caldirola's  expansion  of 
the  distribution  function  is  similar  to  Ginzburg's'^  'treatment,  where  an 
"effective"  energy  loss  for  electrons  by  collisions  is  written 


_eff.  .  R(v) 

®  -  Vw 


Caldlrola  takes  into  account  ionizing  and  recombination  collisions  under  the 
assumption  that  steady  state  conditions  pertain  (which  is  accurate  for  the  excita¬ 
tion  of  an  alrglow,  but  inapplicable  for  ^  breakdown); 


I 


d^  V  =  0 


Caldlrola  has  given  the  following  expression  for  the  inelastic  collision  integral  repre¬ 
senting  ionizing  collisions:  • 


V  V 

=  -Vj^(v)  f(v,c)  +  tj  Vj^(vj^)  f(Vj^,t)  +  n  f(vj^.t) 


where  Vj^(v)  “ 


N.  =  number  of  molecules  per  cc 
A 

.1/2 
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e.  ionisation  energy 


Tj  and  ij  rjl~r^  are  the  parameters  deterciing  the  energy  distribution  between  the 
primary  sesetering  and  secondary  electrons. 
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Carlecon  and  have  nunerically  conpuCed  solutions  o£  the  Boltzmann 

equation  for  the  electron  energy  distrlbatian  in  weakly  ionized  air.  They  have  nade 
the  following  assunptions: 

(1)  That  the  partially  ionized  plastna  is  situated  in  a  static  magnetic  field. 

(2)  That  an  alternating  electric  field  exists  with  electric  vector  perpendicular 
to  the  magnetic  field. 

(3)  That  the  gas  is  homogeneous  and  the  fields  uniform  in  space.  „ 

(4)  The  degree  of  ionization  is  weak  enough  so  that  electron-eleptron  and 
electron-ion  collisions  are  negligible. 

(5)  The  fields  are  such  that  the  average  electron  energy  is  much  larger  than  the 
thermal  energies  of  the  gas  molecules  (and  that  heating  of  the  gas  by  the 
electrons  is  negligible). 

(6)  That  ionization  and  recombination  collisions  are  negligible. 

(7)  That  electrons  may  lose  energy  in  elastic  collisions  and  may  excite  rotational 
vibrational,  and  electronic  degrees  of  freedom. 

Carleton  and  Meglll^  describe  the  various  elastic  and  inelastic  energy 
loss  processes  by  employing  the  known  experimental  cross  sections.  The  gas  was  con¬ 
sidered  to  be  composed  of  N^,  0^  and  0  in  proportions  corresponding  to  an  altitude  of 
about  100  km.  The  total  electron  velocity  distribution  function  was  expanded  in  the 


f(v,t)  =  fQ(v,c)  +  (V'E)  f^(v,c)  +  [v  •  (B  X  E)]  f^Cv.t) 


where  and  f^  provide  bulges  on  the  distribution  in  the  directions  parallel  to  the 
electric  field  and  to  the  Hall  drift,  respectively.  Their  final  result  for  the 
electron  distribution  function  exhibits  a  small  hump  on  the  calculated  curve  near  the 
origin,  which  •..■as  due  to  the  error  introduced  in  die  approximation  wh.ich  treats  the 


O 


rotational  cross  section  as  a  step  function  at  0.02  ev.  It  is  clained  that  the 
calculated  curve  is  inaccurate  beyond  4  or  5  ev  because  higher  lying  states  have 
not  been  included.  If  such  higher  Lying  states  had  been  included,  the  calculated 
curve  would  have  remained  below  the  Druyvesteyn  distribution,  instead  of  rising 
above  it.  It  is  definitively  demonstrated  that  the  resultant  distribution  differs 
greatly  from  the  Maxwellian  over  the  electron  energy  range  2  to  9  ev. 
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Of  the  several  experiments  on  rf  breakdown  currently  being  per¬ 
formed,  that  employing  an  external  rf  source  and  a  transparent  spherical 
cavity  most  nearly  approximates  the  geometry  of  the  real  case  of  break¬ 
down  in  the  atmosphere  produced  by  a  focussed  ground  antennae.  Moreover, 
because  of  the  lack  of  mode  excitation  (and  the  transparent  enclosure), 
it  is  relatively  simple  to  extend  the  scope  of  these  experiments  to 
Include  impurities  and  even  solar  radiation  (according  to  our  information 
on  the  solar  flux  at  the  various  altitudes).  For  cavities  large  enough 
to  avoid  or  inhibit  the  contributions  of  wall  effects,  the  simulated  solar 
flux  will  provide  atomic  species  which  are  impossible  to  reproduce  in 
most  cavity  work,  but  which  are  a  fixture  of  the  upper  atmosphere. 

It  is  to  be  emphasized  that  the  results  of  the  remote  antenna 
breakdown  experiments  differ  from  those  of  comparable  experiments 
employing  microwave  cavities.  However,  circumstances  have  dictated  that 
no  intensive  effort  be  made  to  date  to  set  up  a  detailed  theory  accounting 
for  the  observations.  As  a  result  of  the  special  advantages  of  this 
particular  configuration  (  a  modulated  conical  geometry),  we  suggest  a 
two-fold  program. 

First,  the  experimental  work  should  be  extended  so  as  to  include  all 
combinations  of  pure  gases,  mixtures  of  gases,  varying  pressure,  and 
simulated  solar  radiation. 

Second,  an  attempt  should  be  made  to  employ  the  Boltzmann  Equation 


to  solve  the  tiire  sequence  c£  events  leading  to  breakdown.  Of  course, 
this  second  part  will  be  difficult,  but  for  the  first  time,  there  will 
be  available  reliable  numbers  from  an  experiment  with  a  siiiq)le  and 
pertinent  geometry. 


9.  A  Direct  Measurenent  of  RF  Breakdovn  in  the  Upper  Atmosphere 

It  has  been  pointed  cut  that  there  are  two”b‘asic  obstacles  to  the 
solution  of  the  breakdown  problem  in  the  unbounded  ambient  atmospheric 
environment.  These  are: 

1.  The  simplified  phenomenological  theory  fails  in  the  low  pressure 
range  even  for  pure  gases. 

2.  The  effects  of  an  altering  composition  (this  parameter  is  not 
well  known  above  balloon  heights),  solar  flux,  etc.  cannot  be  reliably 
included  in  either  the  simplified  theory  or  to  date  in  the  more  general 
Boltzmann  approach.  Consequently,  it  seems  that  a  direct  measurement 
upon  the  rf  power  required  for  breakdotm  be  made  in  the  atmosphere, say, 
in  a  parallel  plate  experiment. 

Two  vehicles,  balloons  and  rockets,  immediately  suggest  themselves, 
and  each  has  its  advantages  and  disadvantages.  The  two  basic  parameters 
of  the  vehicle  are  lifting  capability  (therefore  altitude)  and  velocity 
with  respect  to  the  ambient  atmosphere.  Present  balloon  capabilities  are 
limited  to  loads  of  1000  lbs  at  100,000  ft  (30  km),  whereas  boosters  are 
available  to  send  10  Lens  to  such  an  altitude.  On  the  other  hand,  large 
relative  velocities  between  the  carrier  and  tlie  environment  may  lead  to 
some  complication  in  the  ir-tarpretation  of  the  experiment.  .Actually,  it 
is  anticipated  that  this  effect  be  small  particularly  for  pulsed  break- 


Time  does  not  permit  an  examination  of  the  load  required  to  pro- 
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vide  rf  fluxes  of  the  order  of  10  watts/m  .  But,  it  is  strongly  suggested 
here  that  a  detailed  and  thorough  feasibility  study  be  made  if  it  is 
considered  that  exact  Information  on  breakdovm  thresholds  in  the  ambient 
atmosphere  is  required. 


10.  RF  Breakdown  in  Perturbed  Atmospheres 

Once  the  problem  of  rf  breakdown  in  the  ambient  upper  atmosphere  is 
solved,  attention  will  be  focussed  upon  the  obvious  sequel,  rf  breakdown 
in  perturbed  atmosphere.  Here,  there  will  be  two  general  types  of 
perturbations  to  be  considered;  namely,  chemical  seeding  and  nuclear 
detonations. 

In  regard  to  chemical  seeding,  there  will  again  be  two  basically 
different  chemicals  to  be  considered.  Particular  chemicals  have  abnormally 
large  (electron)  attachment  cross  sections  and  therefore  tend  to  inhibit 
breakdown.  The  halogens,  hexafluorides,  etc.  belong  to  this  class  of 
substances.  There  will  be  a  second  class  which  can  lower  the  breakdown 
threshold  because  of  inherently  lower  ionization  potentials.  The  alkali 
metals  are  good  examples  of  these  electron-givers. 

The  nuclear  detonation  -  particularly  at  high  altitudes  -  produces 
a  radical  perturbation  of  the  environment  over  extensive  volumes  of 
the  upper  atmosphere.  Tl'ie  basic  difference  betvreen  the  chemical  seeding 
and  the  nuclear  detonation  is  that  in  the  latter  case  the  perturbation 
is  produced  by  radiation  rather  than  an  injected  impurity.  Nevertheless, 
the  detonation  in  addition  to  the  dominant  effect  of  an  increased  electron 
density  (and  therefore  a  lower  breakdown  field)  can  result  in  an  alteration 
of  Che  composition  through  such  processes  as  dissociative  recoii) inat ion 
and  subseque.nc  atem-ion  or  atom-atom  reactions. 


No  attempt  is  made  here  to  more  than  delineate  the  general  problem 
and  suggest  that  some  consideration  be  given  to  this  area  in  the  future. 


